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Abstract 


In this thesis, we discuss the idea of multigravity, namely the possibility that a signih- 
cant component of gravity that we feel at intermediate distances (1mm < r < 10^®cm) 
is due to massive but ultralight gravitons. We demonstrate how this phenomenon can 
be realized in brane-world models in a spacetime with more than four dimensions and 
discuss how modihcations of gravity at cosmological scales emerge as a consequence. 
Firstly, we consider hve dimensional multigravity models with flat branes. We see how 
the existence of freely moving negative tension branes gives rise to ultralight graviton 
Kaluza-Klein states. Secondly, we study the moduli corresponding to the position of 
the branes and the size of the extra dimension, the radions and the dilaton respectively. 
We show that the radions corresponding to negative tension branes have wrong sign 

kinetic term. We also derive a stabilization condition for the dilaton in a brane model 
with general bulk scalar held dynamics. Thirdly, we show how we can reproduce the 

hve dimensional multigravity models with Anti-de Sitter branes of positive tension 
only. Moreover, we demonstrate how the van Dam-Veltman-Zakharov discontinuity 
can be circumvented at tree level in Anti-de Sitter space. Finally, we discuss how 
the above multigravity models can be replicated in six dimensions with hat positive 
tension branes. This last possibility is the hrst theoretically and phenomenologically 
viable realization of the multigravity scenario. 
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To my parents, Eleni and Yiannis. 
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Chapter 1 
Introduction 


There has been a lot of attention over the past few years in theories where the Standard 
Model (SM) helds are conhned on a three dimensional dynamical hypersurface, a 
three-brane, in a higher than four dimensional spacetime. This scenario is widely 
known as the brane-world conjecture. This idea dates back to the early eighties when 
independently Akama |j^ and Rubakov and Shaposhnikov suggested models of 
our Universe as a topological defect and later Visser and Squires ^ described how 
particles can be gravitationally trapped on the brane. This conjecture has been subject 
of renewed interest in the recent years with the realization that such structures are 
common in string theory where the presence of D-branes [Q where open stings end is 
ubiquitous. 

We can generally divide the recent brane-world models into two classes. The hrst 
one pioneered by Antoniadis, Arkani-Hamed, Dimopoulos and Dvali assumes 

that we live on a 3-brane in a (3 -|- 1 -|- n)-dimensional spacetime with factorizable 
geometry along the extra dimensions. The higher dimensional Planck scale M is then 
related to the 4D Planck scale by Mpj = where W is the compactihcation 

volume. Taking the size of the new n dimensions to be sufficiently large and identifying 
the (4 -|- n)-Planck scale with the electroweak scale (TeV), a hierarchy between the 
electroweak and the Planck scale is generated. This gives an alternative to technicolour 
or low energy supersymmetry resolution to the hierarchy problem. Experimental and 
astrophysical constraints demand that n > 2, M ^ 30TeV and allow new dimensions 
even of submillimeter size. However, a new hierarchy must now be explained, namely 
the ratio of the large compactihcation radius to the electroweak scale. Even for six 
extra dimensions this must be greater than 10^. 

The second class of models which attracted equally considerable attention were the 
ones involving non-factorizable geometry along the extra dimensions and have quite 
distinct phenomenology from the factorizable geometry case. These models come 
under the name of warped brane constructions and are motivated by heterotic M- 
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theory 0.0. ng and its five dimensional reduction In the context of these 

models, one can localize gravity on the brane world having four dimensional gravity 
even with one extra dimension of infinite extent [^, |^ (RS2), or can generate an 
exponential scale hierarchy with a compact two brane model of Planck scale volume 
as it was done in the Randall-Sundrum (RSI) model |^, providing a novel geometrical 
resolution of the hierarchy problem. The scale hierarchy is generated by an exponential 
with respect to the extra dimension warp factor multiplying the Minkowski metric. 
We will briefly review in more details this model at the end of this introduction since 
we will concentrate on warped brane-world models. 

In this thesis we will discuss another interesting possibility which arises from warped 
brane-world constructions, the scenario of multigravity. The proposal of multigravity 
is that the gravitational interaction that we feel in our Universe is the result of the 
combination of massless and massive gravitons or even of massive gravitons only. The 
extra dimensional models give us an infinite tower of massive Kaluza-Klein (KK) 
gravitons after the process of dimensional reduction. However, one should keep in 
mind that Newton’s inverse square gravitational law is very well tested for distances 
that range from the submillimeter regime up to at most lO^^cm that corresponds 
to 1% of the size of the observable Universe (see Fig. |l.l|) . Therefore, the massive 
KK gravitons that will give rise to multigravity should have Compton wavelengths 
bigger than lO^^cm. Additionally, the remaining of the KK tower should either have 
wavelengths entirely below the submillimeter scale or if they span the whole region they 
should have very suppressed couplings. It is more than obvious that these requirements 
are non trivial and cannot be realized in simple KK models since the equidistant 
spectrum that arises in the typical case will be in contradiction with observations. The 
aim of this thesis is to demonstrate how one could take such patterns in warped brane- 
world models. The interesting feature of these multigravity models is that gravity is 
modified at large scales where one starts to be sensitive to the graviton’s mass. We 
should note at this point that this is the first time that modification of gravity at 
cosmological scales is triggered by microphysics considerations which have to do with 
the structure of the extra dimensional space of usually Planck scale volume. These 
scenarios are in principle testable by the Cosmic Microwave Background (CMB) power 
spectrum EH, gravitational lensing [|T^ and gravitational wave observations. 

We will firstly present multigravity models which involve only fiat branes in five 

dimensions. We will begin our presentation by the "H- h" bigravity model [M| where, 

apart from the graviton zero mode, an anomalously light graviton KK state is present. 
This gives the possibility that gravity at intermediate scales (1mm < r < lO^^cm) is 
mediated by both the massless graviton and the ultralight state, leading to the termi¬ 
nology “bigravity”. At large scales were the Yukawa suppression of the massive state 
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Figure 1.1: Our current limit of knowledge of gravity. 


turns on, the gravitational constant decreases to a value which for example in the sym¬ 
metric conhguration is half of its value at intermediate scales. Similar constructions 
of freely moving negative tension branes can also give qualitatively different modi- 
hcations of gravity at large scales. We will briefly review the quasi-localized gravity 
model of Gregory, Rubakov and Sibiryakov (GRS) which has no normalizable zero 
mode but gravity at intermediate distances is generated by a “resonance”-like [P5| , |2^ 
coupling of the lowest part of the continuum of the KK spectrum. In this construction, 
the gravitational law at large scales changes from four dimensional to hve dimensional. 
Although the bigravity and the quasi-localized gravity models look quite different, they 

. We will explicitly show 


are connected through a more general ” -\ - \-" model |25 


how the spectrum and the gravitational law interpolate between the ones obtained in 
the above cases. The change of the nature of the gravitational law is observed also in 
the the crystalline universe model [^] which we will discuss afterwards. In this case 
gravity at intermediate distances is generated by the hrst low lying band of KK states 
and at large scales gives again hve dimensional law of gravity. For an alternative way 
of modifying gravity at ultralarge distances see [^,1^. At this point we will present 


a potential danger arising from the fact that massive gravitons have extra polariza¬ 
tion states in comparison with the massless one. This is clearly seen by the difference 
in the tensor structure of the massive and massless spin-2 propagators in hat space 
(van Dan-Veltman-Zakharov discontinuity) 


i|,|30|. Observationally this implies that 
the massive graviton proposal disagrees with well tested predictions of Einsteinian 
gravity [Q . 


So far our discussion has considered the gravitational force due to tensorial gravi¬ 
tational modes. However, in the gravity sector there exist moduli corresponding to the 
huctuation(s) of the brane(s) or the size of the extra dimension, the radion(s) [p^[33[ 
and the dilatonQ respectively, which change the phenomenology of these mod¬ 

els. This will be the subject of chapter 3 where we will derive the ehective action de¬ 
scribing these excitations . For this purpose we will hrstly consider a general three 
three-brane model with arbitrary bulk cosmological constants between the branes. We 


^This is not to be confused with the dilaton field that is present in the spectrum of closed strings. 
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will calculate the kinetic energy for the moduli of the system and see that the freely 
moving (i.e. not on an orbifold hxed point) negative tension branes have wrong sign 

This plays an important role in the context of multigravity 


kinetic term 31 


models where this ghost radion(s) cancel the troublesome extra polarization states of 
the massive gravitons restoring the phenomenologically favoured tensorial structure of 
the graviton propagator At large distances the ghost radion(s) dominate 

giving scalar antigravity 0 in all the above constructions instead of a reduction of 
the Newton’s constant or a change of the nature of the gravitational law. However, 
the very fact that there are physical ghost(s) in the spectrum makes these models 
unacceptable. Note that the essential characteristic of all these models is the bounce 
of the warp factor generated by the presence of the moving negative tension branes. 
This association between the bounce in the warp factor and presence of negative ten¬ 
sion branes is an unavoidable feature in flat (five dimensional) brane models and is 
possibly linked to the fact that the weaker energy condition is violated at the position 
of the bounce In the remaining of the chapter we will discuss the dilaton 


for backgrounds with identical four-geometries and with bulk dynamics involving ar¬ 
bitrary scalar helds minimally coupled to gravity [0. In order to do this we will 


briefly review the brane-world sum rules derived by Gibbons, Kallosh and Linde [41 


Then we will find a generic formula relating the dilaton mass with the sum of the 
tensions of the branes and the curvature of the four-geometries, from which we can 
see immediately the necessary condition for stabilizing the overall size of the system. 

The flat multigravity models as we have already said suffer from the presence of 
ghost scalar mode(s) which renders the theory unacceptable. In chapter 4 we will 
present a way out of this difficulty by considering AdS 4 branes instead of flat branes. 
In this case the warp factor has a " cosh^' form which naturally generates a bounce 
without any need of floating negative tension branes. Indeed, it is straightforward 
to replicate the bigravity and the crystalline model in this framework. The bigravity 
model is converted to a " +" model and the crystal model to an infinite array of 

branes. By sending the second brane in the ”++” model to infinity, one obtains the 


locally localized model of where gravity is mediated only by the ultralight 

state since the zero mode becomes non-normalizable. It is interesting to note how these 
models confront the relativistic predictions of Einsteinian gravity since there are no 
ghost radion(s) to restore the graviton propagator structure to the phenomenologically 
favoured one. We will show that in {A)dS space the van Dam-Veltman-Zakharov no- 
go theorem can be evaded at tree level for the case that m/H —> 0, where m is 
the mass of the graviton and H the “Hubble” of the {A)dS space, as it was shown 
^,T7,^,4^. The extra polarizations of the massive gravitons in this case are 


m 
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very weakly coupled and decouple in the massless flat limit. A precise realization 
of this case is the " + +" bigravity model. The AdS brane models, even though 
they solve the theoretical difficulty of the ghost radion(s) of the flat brane models, 
face phenomenological difficulties. The presence of a remnant negative cosmological 
constant is at odds with observations, and furthermore it turns out that all large 
scale modifications of gravity predicted by these models, are hidden behind the AdS 
horizon. Furthermore, as we will discuss at the end of the chapter an old result by 
Vainshtein suggests that the discontinuity of the graviton propagator in flat space 
is not directly observable by the bending of the light by the Sun as it was claimed 
in because non-perturbative effects become important and make the massless 

limit continuous. 

In the light of this encouraging result it seems appropriate to look for multigravity 
models of flat branes without ghost fields. As we have discussed, in five dimensions it 
is impossible to have flat brane multigravity models without negative tension branes. 
For this reason, in chapter 5 we will consider models in six dimensions |^. We 


explicitly show that the difficulties of the five dimensional models can be evaded in 
six dimensions. It is possible to construct multigravity models with only flat positive 
tension branes. The branes which localize gravity in this setup are four-branes, but 
one of their dimensions is compact, unwarped and of Planck length. Thus, in the low 
energy limit the spacetime on the brane appears four dimensional. In order that these 
constructions are realized it is crucial that the tensions and/or the bulk cosmological 
constant are anisotropic. The five dimensional constructions in this setup come into 
two types. One of them involves conical singularities at finite proper distance from the 
four branes. These conical singularities support three-branes which can be of positive 
tension if one has an angle deficit, of zero tension if one has no angle deficit and of 
negative tension if one has an angle excess. The other type of the multigravity models 
has no conical singularity at all. It is crucial to note that there is no five dimensional 
effective theory for these constructions, otherwise one would get all the problems faced 
in the five dimensional constructions. The low energy effective theory is directly four 
dimensional. These six dimensional models give for the first time a theoretically and 
phenomenologically viable realization of multigravity. 

In the final chapter, we summarize the multigravity models that we have presented 
in this thesis and discuss questions for future investigation. 

Throughout this thesis we use the mostly plus metric convention with signature 
being and we use the following convention for the Riemann tensor 

R\mn = ^m'^an - Sn'^am + ^an'^mh - ^am'^nh where the Christoffel symbols are 
^MN = ^9^^ {^MgAN + dj^gAM “ We work in the units where h = c = 1. The 

four dimensional Planck mass and Newton’s constant are related as Gj/ = 327rMp[. 
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Figure 1.2: The Randall-Sundrum model with one and one " brane at the hxed 
points of a jZ^ orbifold. 

1 The Randall-Sundrum model 

Since, throughout this thesis we will study warped models it is instructive to briefly 
present the Randall-Sundrum model. This construction consists of two parallel three- 
branes sitting on the hxed points of an /Z 2 orbifold in hve dimensional spacetime 
which is essentially a slice of AdS^ (see Fig. p..2|) . The negative bulk cosmological 
constant and the tensions of the branes are hne tuned in such a way that the four¬ 
dimensional slices of the hve-dimensional spacetime, and thus the two branes, are hat. 
The metric that describes this conhguration has the characteristic warped form: 

ds^ = e~'^^^'^^r]^ydx^dx’^ + d'lf ( 1 . 1 ) 

If one now considers a scalar held, in particular the SM Higgs, conhned on one of 
the two branes, then the minimal coupling to the metric will result to an action: 


d X ^ y/ Qind ^ \ 

- si 

e-*'^\yi\ 1 

^ J 


( 1 . 2 ) 


where is the induced metric on the brane at position yi. After a held rescaling 
H —>■ which brings the kinetic term to be canonical, the action reads: 

^H = J (1.3) 

which clearly shows that the vacuum expectation value (vev) of the Higgs is given as: 

-k\yi\ 


V = VqC 


(1.4) 






1 The Randall-Sundrum model 
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m 


Figure 1.3: The graviton (blue line), first (red line) and second (red dashed line) KK 
states wavefunctions in the RSI model. On the right the pattern of the spectrum is 
sketched. 


Assuming that all hve dimensional parameters and thus uq are of the order of the 
Planck scale, then a Higgs on the first brane will have a vev of the order of the Planck 
scale, whereas a Higgs on the second brane will have exponentially reduced vev. It was 
exactly this generation of exponentially different mass scales that was utilized in the 
Randall-Sundrum scenario to resolve the Planck hierarchy problem by placing the SM 
on the second brane (TeV brane). A compactification radius of only some 35 times 
larger than the Planck length was sufficient to bring the Higgs vev to the TeV scale. 

The phenomenology of this model has been extensively explored over the past years 
(see for a comprehensive account). In the TeV hierarchy scenario, the KK tower 
of the spin-2 graviton resonances starts from the TeV scale with approximately TeV 
spacing between them giving rise to characteristic signals in high energy colliders. The 
graviton and KK states wavefunctions {i.e. their profiles along the extra dimension) 
are shown in Fig. |1.3| . The massless zero mode wavefunction follows the warp factor and 
is localized on the positive tension brane, while the wavefunctions of the remaining 
tower are suppressed on the positive tension brane and enhanced on the negative 
tension one. We will return to this picture later and compare it with the considerably 
different pattern that we get in the multigravity models. 










Chapter 2 


Flat brane multigravity models in 
five dimensions 


In this chapter we will see how we can realize the multigravity scenario, as it has been 
described in the introduction, with flat multi-brane conflgurations in flve dimensions. 
We will see that these models are infested with the presence of freely moving negative 
tension branes. We will discuss the problems arising from the negative tension branes 
in the following chapter. Here, we will concentrate on the gravitational potential 
coming from the tensor fluctuations of the background metric and discuss the general 
characteristics of the multigravity models. Although flat brane multigravity scenarios 
have difficulties, they are very easy to handle toy models and can be used to understand 
physics of more complicated but realistic multigravity models. 

In the following we will discuss the " H-bigravity modes, the GRS model of 

quasi-localized gravity and the relation between them as well as the crystal multigrav¬ 
ity model. 


1 The " + —h bigravity model 


We begin our presentation of flat multigravity models in flve dimensions with the 

n 


H- V" bigravity model |^|. This model is a straightforward extension of the RSI 

model with one more positive tension brane. In more detail, we have three parallel 
three-branes in an AdS^ space with cosmological constant A < 0. The fifth dimension 
has the geometry of an /Z 2 orbifold and the branes are located at Lq = 0, Li and 
L 2 where Lq and L 2 are the orbifold fixed points (see Fig. |2.1| ). This model has three 
parameters, namely Li, L 2 and A. For our present purposes we will concentrate in 
the symmetric configuration where the intermediate brane lies in the middle of the 
orbifold, i.e. L 2 = 2Li = 21, leaving two parameters /, A. Firstly, we consider the 
branes having no matter on them in order to And a suitable vacuum solution. The 
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Figure 2.1; The " H-bigravity model with two branes at the hxed points of 

a /Z 2 orbifold and a brane in between. 


action of this setup is: 


Sg 


o2l 


d^x / dy^/^{2M^R- K] 


'-21 



( 2 . 1 ) 


where is the induced metric on the branes, Vi their tensions and M the hve 
dimensional fundamental scale. The integration is performed throughout the orbifold, 
so there is no need to introduce Gibbons-Hawking boundary terms. The Einstein 
equations that arise from this action are: 

Rmn - \gmnR ("aGma- + - U) j (2.2) 

At this point we demand that our metric respects four dimensional Poincare in¬ 
variance. The metric ansatz with this property is the following: 

ds^ = -|- dy^ (2.3) 

Here the warp function <j{y) is essentially a conformal factor that rescales the 
four dimensional component of the metric. A straightforward calculation gives ns the 
following differential equations for (j{y)\ 


{a'f = e 

(2.4) 


(2.6) 
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where k = 


-A 

24M3 


is a measure of the curvature of the bulk. 


The solution of these equations consistent with the orbifold geometry is precisely: 


a(p) ^t{l- ||j/| - i|} 


( 2 , 6 ) 


with the additional requirement that Vq = V 2 = —Vi = V. Furthermore, four di¬ 
mensional Poincare invariance like in the RSI model requires the hne tuning relation: 


V = 


A 

k 


(2.7) 


If we consider massless fluctuations of this vacuum metric as in 16 and then 


integrate over the hfth dimension, we hnd the four dimensional fundamental scaleQ 
M*, the hve dimensional fundamental scale M and the curvature of the bulk k are 
related as: 


2M3 


= — [1 - e-“‘] 


( 2 , 8 ) 


We do not identify the mass scale M* with the four dimensional Planck scale for 
reasons which will be obvious in the following. The above formula tells us that for large 
enough kl the three mass scales M*, M, k can be taken to be of the same order. Thus 
we take k ~ 0{M) in order not to introduce a hierarchy, with the additional restriction 
k < M so that the bulk curvature is small compared to the hve dimensional Planck 
scale and we can trust our solution. Furthermore, if we put matter on the branes all 
the physical masses m on the brane will be related to the fundamental mass parameters 
mo ~ M of the hve dimensional theory by the conformal warp factor: 


m = e 


(2.9) 


Thus, in the symmetric conhguration which we are examining, the mass scales 

on the hrst and third brane will be of the order of the fundamental scale M while 
on the intermediate negative tension brane they will be exponentially suppressed. It 

is obvious that in an asymmetric conhguration we could obtain an electroweak scale 
hierarchy a la Randall-Sundrum even on a positive tension brane. For this case we 
would only need 2Li — L 2 ~ in order to get ~ O (see pl|,[53[]). 

To determine the phenomenology of the model we need to know the KK spectrum 
that follows from the dimensional reduction. This is determined by considering the 
linear huctuations of the metric of the form: 


^ [e y)] dx^dx'' + 


( 2 . 10 ) 


^The factor M* multiplies the four dimensional Ricci scalar in the zero mode truncated Lagrangian. 
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Here we have ignored the dilaton and radion modes which will be the snbject of 
the following chapter. We expand the held h^u{x,y) in graviton and KK states plane 
waves according to: 

CX) 

where {d^d^ — ml) = 0 with = 0. The latter conditions for the 

zero mode correspond to a choice of gauge hxing. The function will obey the 

following second order differential equation: 

+ 2 V - 2 k life) + Hy - 21 ) - ife - 0 - S(y + /)]U'">fe) = 

( 2 . 12 ) 

After a change of variables and a redehnition of the wavefunction the above equation 
reduces to an ordinary Schrodinger-type equation: 

j-isj + r(fe| (2,13) 

with potential: 

= or / 312 “ -<5(^ + ^z)] (2.14) 

2g{z) 

The new variables and wavefunction in the above equation are dehned as: 


z = sgn{y) 


2gfci _Q'2kl — k\y \ _ ^ 

k 

e^\y\ — l 
k 


, / < bl < 2/ 
Ay\<i 


(2.15) 


4zW(^) = ^W(^)e-/2 (2.16) 

and the function g{z) as g{z) = k{zi — \ \z\ — z)} + 1, where zi = z{l). The change of 

variables has been chosen so that there are no hrst derivative terms in the differential 
equation. Furthermore, in this coordinate system it can be easily seen that the vacuum 
metric takes the conformally flat form: 

ds^ = {rjavdx^dx^ + dz^) (2-17) 

g[zy 

This is a quantum mechanical problem with 5-function potentials and an extra 1/g'^ 
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Figure 2.2: The potential V(z) in the " H- h" bigravity model. The (5-function wells 

correspond to branes, while the (5-function barriers to branes. 


smoothing term (due to the AdS geometry) that gives the potential a double “volcano” 
form (see Fig. |2.2|) . An interesting characteristic of this Schrodinger equation is that it 
has a “supersymmetric” form which guarantees that it has no tachyonic modes: 


with A=^--^ 
^ ^ n V / dz 2 g 


(2.18) 


Since the system is compact, this differential equation always gives rise to a massless 
zero mode which reflects the fact that Poincare invariance is preserved in the four 
dimensional spacetime. Its wavefunction is given by: 


4 ^( 0 ) = 


A 


W)fi^ 


(2.19) 


The normalization factor A is determined by the requirement 



^( 0 )( 2 ) 


1 , 


chosen so that we get the standard form of the Fierz-Pauli Lagrangian. 

In the symmetric case that we are examining, the potential, and thus the zero 
mode wavefunction, is symmetric with respect to the second brane. This immediately 
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implies that the next level will have a wavefunction antisymmetric with respect to 
the second brane. This is exactly what happens. The solution for the graviton KK 
modes is given in terms of Bessel functions. For y lying in the regions A = [0, /] and 
B = [/, 2/], we have: 


\|/W 


A 

B 



( 2 . 20 ) 


The boundary conditions (one for the continuity of the wavefunction at zi and three 
for the discontinuity of its first derivative at 0, zi, Z 2 ) result in a 4 x 4 homogeneous 
linear system which, in order to have a non-trivial solution, should have vanishing 
determinant: 


"'i(f) 

0 

my) 

■h my 

Here we have suppressed the subscript n on the masses m„. This is essentially the mass 
quantization condition which gives the spectrum of the KK states. For each mass we 
can then determine the wavefunction with the same normalization as the zero mode. 

The exact calculation of the wavefunction of the hrst graviton KK state verihes our 

previously stated intuition. Another way to understand the form of the wavefunction is 
that since we have two delta functions in the potential each of which supports a bound 
state, there should be a massive KK state that resembles the zero mode. Indeed, if one 
calculates the wavefunctions of the higher levels it is evident that the hrst KK state 
has a special characteristic. It is the only one that is peaked in the positive tension 
branes, whereas all the higher levels are suppressed on the positive tension branes and 
are spread in the bulk (see Fig. p.3|) . This is in contrast to the RSI case where the 
wavefunctions of all the KK states are spread in the bulk (see Fig. |1.3|) . Moreover, 
since the difference of the modulus of the wavefunctions of the hrst special KK state 
and the zero mode comes from a very small region of the r/-space around the position 
of the minimum of the warp factor, we expect its mass to be anomalously light in 
comparison to the masses of the remaining of the KK tower. 

This is a generic characteristic of potentials which have the shape of a bounce. 
In hve dimensions with hat branes the only possibility of obtaining a bounce is by 
introducing negative tension branes. This is because the warp factors generated by 
hat positive tension branes always have the form of a falling exponential. The only way 


K. (f) 0 0 

0 Ji (¥) n (t) 

n my -h my v my 
nmy -^2my -y^my 
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Figure 2.3: The graviton (blue line), first (red line) and second (green line) KK states 

wavefunctions in the symmetric " H- h" model. The wavefunctions of the zero mode 

and the second KK state have discontinuities on the branes. On the right the 
pattern of the spectrum is sketched. 

of pasting two warp factors which are falling exponentials is through a discontinuity 
generated by a negative tension brane. 

If we dehne x = kl, as the distance between the negative and positive branes 
measured in units of k, then for x S> 1 we may obtain a reliable approximation to its 
mass by using the hrst terms of the Bessel power series. Then we hnd that the hrst 
KK mode has mass given by: 

mi = 2V2ke-‘^^ (2.22) 

The masses of the other KK states are found to depend in a different way on the 

parameter x. Analytically we hnd that for x 3> 1 the masses are: 

mn+i = inke~'" , n = l,2,... (2.23) 

where ^ 2 i+i is the (i + l)-th root of Ji(x) {i = 0,1, 2,...) and ^ 2 i is the i-th root of 
J 2 (x) {i = 1,2,...). Using the asymptotic values of the Bessel roots, we hnd that the 
level splitting for high n levels approaches a constant value: 

Am = (2.24) 

The coupling also of the hrst KK state to matter on the branes is very diherent 
from the one of the remaining tower. To see this, we consider a minimally coupled 
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to gravity matter Lagrangian Cmi^j, where we denote with the matter fields 
living on the branes. The total action of the system is the sum of the gravity action 
(|2.1|) plus the matter action; 


S = SG + Y,j d^xdy^f^)CU^j, 


(2.25) 


The effective four dimensional Lagrangian taking into account only the graviton 
perturbations, has the following form: 


C = Cmm,r]) + 2M- 




26) 


where Cpp is the Pauli-Fierz Lagrangian for the n-level graviton mode. The matter 
fields have been rescaled to make them canonically normalized and the 

energy momentum tensor Tj^'^ for matter localized on the i-th brane is defined with 
respect to the rescaled fields <h^. By construction the induced background metric on 
the branes is the fiat four dimensional Minkowski metric Finally, defining the 

canonically normalized graviton fields: 




(2.27) 


the Lagrangian reads: 




2vW 


(2.28) 


Thus, the gravitons have the following couplings to matter on the positive tension 
branes: 

{z+) g 


^T). - 


2V2IP 


(2.29) 


In more detail we find that the zero mode and the first ultralight KK state couple 
with the same strength on the positive tension branes: 

1 


2V2M, 


Go — 0,1 — 


(2.30) 
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while the coupling of the remaining tower is exponentially suppressed: 


®n+l 


1 e-^ 

2V2M, VTEyTTEy 


n = 1, 2, 3,... 


(2.31) 


The reason for the strong coupling of the hrst KK state in comparison with the 
remaining tower is readily understood because, being dominantly a bound state of the 
volcano potential on the positive tension brane, it is largely localized on it. 

Let us now demonstrate how the bigravity scenario can be realized in this model. 
We have already seen from (|2.22|) and (|2.23| ) that, for large x, the lightest KK mode 
splits off from the remaining tower. This allows the possibility that the Newtonian 
gravity that we experience at intermediate distances (1mm < r < 10^®cm), is due to 
both the massless graviton and the ultralight KK state. The gravitational potential is 
computed by the tree level exchange diagrams of the massless graviton and KK states 
which in the Newtonian limit is: 


V(r) 


E 


n=0 


al e 


dvr r 


(2.32) 


where a„ is the coupling (|2.29|) and n = 0 accounts for the massless graviton. The 
summation stops at some very high level N\ with mass of the order of the cutoff scale 
~ M. In the bigravity scenario, at intermediate distances ^ r -C the hrst 
KK state and the massless graviton contribute equally to the gravitational potential: 


Vid{r) 


1 n 1 \ 

327iM^ \r ry 


r 


(2.33) 


where Gn = joyyp {Mpi = ^). For long distances r > the Yukawa suppression 
effectively reduces gravity to half its strength: 


Vidir) 


1 (I 0\ 

327rM^ \r ry 


Gn 1 

2 r 


(2.34) 


Astronomical constraints and the requirement of the observability of this effect 
demand that for k ~ Mpi we should have x in the region 65 — 70. The observational and 
experimental constrains in the context of the graviton KK modes are m < 10“^^eV or 
m > 10“'^eV for modes with Planck suppressed coupling. If we choose rui ~ 10“^^eV, 
from (|2.22|) and (|2.23|) we hnd m 2 ~ 10“^eV which comfortably satishes the Cavendish 
bound (see Fig.|^). 
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0 



r 


Figure 2.4: Exclusion regions for the bigravity case and correlation of the hrst two KK 
states. 


However, the Cavendish bound in this case is not even applicable because the 
coupling of the remaining KK tower is exponentially smaller than the usual Planck 
suppressed coupling. To see this let us calculate the correction to the gravitational 
law at distances r < The Yukawa interactions of the remaining KK states then 

turn on and will give rise to a short distance correction. This can be evaluated by 


using the asymptotic expression of the Bessel functions in ( 2.31 ) since we are dealing 
with large x and summing over a very dense spectrum, giving: 


AK,,(r) 


Gn 

~T 


E 


n=2 


kn 

2e^ 



^-rrinr 


r 


(2.35) 


At this point we exploit the fact that the spectrum is nearly continuum above m 2 
and turn the sum to an integral with the hrst factor in ( p.35| ) being the integration 
measure, i.e. = '^Am —J dm from (|2.24|) . Moreover, we can extend the 

integration to inhnity because, due to the exponential suppression of the integrand, 
the integral saturates very quickly and thus the integration over the region of very 
large masses is irrelevant. The resulting potential is now: 


A14Yr) = 


Gn , m e~ 
k 2k r 


(2.36) 


'm 2 


The integration is easily performed and gives: 


AKrf(r) = 


Gn 1 + rn2r _ 


m2V 


2r {kry 


(2.37) 


We see these short distance corrections are signihcant only at Planck scale distances 
~ k~^. Thus, the lower bound of the Cavendish experiments is not important in 
the symmetric conhguration simply because although the low lying KK states have 
Compton wavelength at the millimeter scale, they are however very weakly coupled. 
Finally, let us comment about the realization of the bigravity scenario in asymmet- 
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Figure 2.5: The GRS model with one "+" brane at the hxed point of an FCIZ 2 orbifold 
and a " — 1/2" brane at a distance x = kL. The curvature of the bulk between the 
"+" and " — 1/2" branes is k and zero in the outside region. 


ric conhgurations where the positive tension branes are not scale equivalent. In this 
case a whole rich of bigravity scenarios emerge. Depending on which positive brane 
we live in, we can have either a dominant massless graviton component or a dominant 
massive graviton component. In the hrst case the correction to Newtonian gravita¬ 
tional law at large scales will be negligible. In the second case, where the massive KK 
state is the dominant source of Newtonian gravity, one would see gravity to essentially 
switch off at large scales. 


2 The GRS model and its relation to the " H- v" 

bigravity model 

An other interesting model which also has the feature of predicting modihcations of 
gravity at ultralarge scales is the one by Gregory, Rubakov and Sibiryakov (GRS) 

In this section we will briefly review the physics of the GRS model and demonstrate 

how it is related to the " -I - h" bigravity model. The GRS model modihes the RS2 

model by adding a negative tension brane of half the tension of the positive tension 
one at some distance L in the bulk. The space between the positive and negative 
tension branes is a slice of AdS^, while in the outside region it is flat (see Fig. p.5|) . 
The background metric for this setup has the four dimensional Poincare invariant form 
(|2.3|) with the warp factor being: 




k\y\ ,\y\<L 
kL ,\y\>L 


(2.38) 


This model does not have a normalizable massless graviton because the volume of 
the system is inhnite, but has a continuum of KK states that starts from zero. Despite 
the absence of a zero mode, it generates four dimensional gravity at intermediate 
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distances due to the “resonance”-like behaviour of the values of the wavefunctions of 
the KK states continuum at the position of the positive tension brane. In that respect, 

gravity in this picture appears to be quasi-localized. The coupling for small masses of 
the KK states has a Breit-Wigner form: 


a{rnf 


A 

m? + (r/ 2)2 


(2.39) 


where T = 4/ce is the width of the resonance. This feature of the coupling gives for 
intermediate distances r -C T”^ normal four dimensional gravitational law: 


V,d{r) 



(2.40) 


where Gn = distances r S> T however, the gravitational law 

becomes hve dimensional: 


Vid{r) 


4Gjv 1 
ttT r2 


( 2 . 41 ) 


If we have L k, then T”^ ~ lO^^cm, so this scenario cannot be excluded by 

current observations. The short distance corrections, as in the " -I- V" model, appear 

only for r ^ k~^ and read: 


AKrf(r) 


Gn 1 

2 r 


(2.42) 


Although the GRS and the "H-bigravity models look quite different, they share 

the key element of freely moving negative tension branes. It was shown in ||^ that 

these two models are the limiting cases of a more general "H-multi-brane model 

that interpolates between them. This model can be readily derived from the bigravity 
model by cutting in half the intermediate negative tension brane, i.e. instead of 
having one negative tension brane one can take two negative tension branes of half the 
tension of the original one, and considering flat space between the two new negative 
tension branes. If the second positive tension brane is taken to inhnity together with 
one of the negative tension branes, we shall get precisely the GRS picture. 

Let us discuss the four-brane " H-model in more detail. It consists of 

a hve dimensional spacetime with orbifold topology with four parallel three-branes 
located at Lq = 0 , Li, L 2 and L 3 , where Lq and L 3 are the orbifold hxed points 
(see Fig. |2.6|) . The bulk cosmological constant A is negative between the positive and 
negative tension branes and zero between the two negative tension branes. The model 
has four parameters, namely Li, L 2 and L 3 and A. For our present purposes we 
consider the symmetric conhguration, leaving three parameters, /, and A where 
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= kC 



- 1/2 - 1/2 

Figure 2.6: The " H- \-” configuration with scale equivalent branes. The 

distance between the and " — 1/2" branes is I = Li = — L 2 while the distance 

between the " — 1/2” branes is F = L 2 — Li. The curvature of the bulk between the 
”+” and ” — 1/2” branes is k and zero between the ” — 1/2” branes. 


I = Li = — L 2 and F = L 2 — Li. In the absence of matter the model is described 

by ( p.lD with 


f 0 ,y e [Li, L 2 
\a ,ye[0,L^][j[L2,L,] 


(2.43) 


By considering the ansatz (^) the warp function a^y) must satisfy: 


a 


n 


[<J 


/\2 


F - A 

i 

f 0 ,y e [Li, L 2 ] 

\ k^ ,ye[0,L,][j[L2,L,] 


(2.44) 

(2.45) 


where k 


-A 

24M3 


is a measure 


of the bulk curvature. The solution for ?/ > 0 is: 




ky ,|/e[ 0 ,Li] 

kl ,ye[Li,L2] 

kl + k{l + F — y) ,y E [L 2 , T 3 ] 


(2.46) 


with the additional requirement that Vq = V 3 = —2Vi = — 2 V 2 = V. Furthermore, 
four dimensional Poincare invariance requires the fine tuning relation: 


V = 


A 

I 


(2.47) 


In order to determine the mass spectrum and the couplings of the KK modes we 
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consider linear fluctuations around the background metric as in ( [1.10| ). Following the 
same procedure we And that the function obeys a Schrodinger-like equation 

(|2.13|) with potential V{z) of the form: 


= 




8[9{z)? 

3k 


[9{z) - 9{z - zi) + 9{z - Z 2 ) - 9{z - zs)] 


2 g{z) ~ ~ ~ ^ ~ 


The conformal coordinates now are given by: 


(2.48) 


z = < 


k 
kl 


{y — /)e^ + 


gKt_l 


, 1/ G [0, Li] 

) 2 / ^ [-^ 1 ) L 2 ] 


_ lg2fci+fci_g ky _|_ _|_ 2gfcZ _ 1 ^ 


(2.49) 


with 


kz -\-l ,ze[ 0 , zi] 

g{z) = { kzi + l ,^€[^ 1 , 2 : 2 ] 

k{z 2 - z) + kzi + 1 ,ze lz 2 , 2 : 3 ] 


(2.50) 


where zi = z{Li), Z 2 = ^(^ 2 ) and Z 3 = z^L^). 

The potential (p.48|) again gives rise to a massless graviton zero mode whose wave- 
function is given by (|2.19|) with the same normalization convention. Note, however, 
that in the limit /_ —> 00 this mode becomes non-normalizable (GRS case). The 
solution of the Schrodinger equation for the massive KK modes is: 


r A ] 


' [ n iygi^)) + AJ 2 {yg{z))] 


> =N^< 

Bi cos{ninZ) + B2 sin(m„2;) 

1 c J 


, sl^-^[C,Yygfg{z))+C,Jgz^g{z))] ^ 


(2.51) 


where A = [0, zi], B = [^ 1 , 2 : 2 ], and C = [z 2 ,Z 3 ]. We observe that the solution in the 

first and third interval has the same form as in the " -I- V" model. The new feature 

is the second region (flat spacetime). The coefficients that appear in the solution are 
determined by imposing the boundary conditions and normalizing the wavefunction. 

The boundary conditions (two for the continuity of the wavefunction at zi, Z 2 and 
four for the discontinuity of its first derivative at 0, Zi, Z 2 and Z 3 ) result in a 6 x 6 
homogeneous linear system which, in order to have a non-trivial solution, should have 
vanishing determinant. It is readily reduced to a 4 x 4 set of equations leading to the 
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quantization condition; 



— cos(m 2 ;i) 

(fs) - ^^ 7 * (fs) 

sin(m 2 ;i) 

0 

— sin(m 2 ; 2 ) 

0 

003(7712:2) 


— sin(mzi) 

— cos(mzi) 

cos(mz2) 

sin(mz2) 


0 

0 

(f5) - ifg) 

-Y2{^g) + ^J2{^9) 


= 0 (2.52) 


with g = g{zi) = g{z 2 ). Here we have suppressed the subscript n on the masses mn- 


2.1 Spectrum Interpolation 

The above quantization condition provides the mass spectrum of the model. It is 
convenient to introduce two dimensionless parameters, x = kl and = kC (see 
Fig. p.6|) and we work from now on with the set of parameters x, X- and k. The mass 

spectrum depends crucially on the distance x_. We must recover the " H- \-" bigravity 

model spectrum in the limit x_ —> 0, and the GRS spectrum in the limit X- —>■ oo. 
From the quantization condition (|2.52|) it is easy to verify these features and show how 

the " H-spectrum smoothly interpolates between the " H- C' bigravity model 

and the GRS one. It turns out that the structure of the spectrum has simple a:_ and 
X dependence in three separate regions of the parameter space: 

The " H-h" Bigravity Region 

For ^ 1 we hnd that the mass spectrum is effectively -independent given by the 
approximate form: 


mi = 2V2ke-^^ (2.53) 

m„+i = ^nke-^ n= 1,2,3,... (2.54) 

where ^ 2 i+i is the {i + l)-th root of Ji{x) {i = 0,1, 2,...) and ^ 2 i is the Tth root of 
J 2 {x) {i = 1, 2, 3,...). As expected the mass spectrum is identical to the one of the 

" H- \-" bigravity model for scale equivalent positive tension branes. The hrst mass 

is manifestly singled out from the rest of the KK tower and for large x leads to the 
possibility of bigravity (see FigCT- 
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Figure 2.7: The behaviour of the mass of the first five KK states in the three regions 
of simple x, dependence. The first dot at zero stands for the graviton. 


The Saturation Region 

For 1 <C x_ <C we find a simple dependence on x_ given by the approximate 
analytic form: 

-2x 

mi = 2k - (2.55) 

^/xT 

^-X 

run+i = mik - n = 1,2,3,... (2.56) 

x_ 

As x_ increases the first mass decreases less rapidly than the other levels (see Fig.|2.7|). 


The GRS Region 

For x_ S> the first mass is no longer special and scales with respect to both x and 
X- in the same way as the remaining tower: 


nin = mrk - n = 1, 2, 3,... 

X- 


(2.57) 


The mass splittings Am tend to zero as x_ —> oo and we obtain the GRS continuum 
of states (see Fig. p.7|) . 
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2.2 Multigravity 

Armed with the details how the spectrum smoothly changes between the " H- h" 

bigravity model (a;_ = 0) and the GRS model (a;_ —*■ cxd), we can now discuss the 
possibilities for modifying gravity at large distances. The couplings of the KK states 
with matter on one of the scale equivalent positive tension branes are readily calculated 
by (PI): 

3/2 

(2.58) 


CLri. 


g(0) 

2V2 I M 


Bigravity Region 

As X- —> 0, we are in the bigravity region which has been extensively explored in the 
previous section. 


The GRS Region 

In the GRS limit, X- 3> we should reproduce the “resonance”-like behaviour of 
the coupling in the GRS model. In the following we shall see that indeed this is the 
case and we will calculate the hrst order correction to the GRS potential for the case 
that x_ is large but hnite. 


For the rest of the section we split the wavefunction ( |2.51| ) in two parts, namely 
the normalization and the unnormalized wavefunction i.e. = 

. The former is as usual chosen so that we get a canonically normalized 

Pauli-Fierz Lagrangian for the four dimensional KK modes (see ( |2.27]) ) and is 
given by: 

1/2 


Nf = 




(2.59) 


Zl 


The value of ip^'^\z) on the left positive tension brane is, for <C /c: 


^h(o) = 






(2.60) 


It is convenient to split the gravitational potential given by the relation ( p.32|) into 
two parts: 


V{r) = - 


8M3 


E 

n=l 


Nk 


-A'>(1)(0) - g/p E - 


(2.61) 


n=Nx 
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where we denote by Nx_ the level after which the coupling of the KK states becomes 
linear to their mass. It corresponds to a mass scale about = ke~^. As we shall 
see this separation is useful because the hrst Nx_ states give rise to the long distance 
gravitational potential Vid while the remaining ones will only contribute to the short 
distance corrections AVsd- 

• Short Distance Corrections 

We hrst consider the second term. The normalization constant in this region is 
computed by considering the asymptotic expansions of the Bessel functions with ar¬ 
gument iAlTlhL^ It is calculated to be: 


= 


TT^m^ 


32k^g{zi)x- 


1 + ^ 
X — 


(2.62) 


If we combine the above normalization constant with the unnormalized wavefunc- 


tion (12.601) , we hnd: 


AWrf(r) 


1 

327rM^ 


E 


n=Nx 


kn 


mn 

2k 


^-rrinr 


r 



(2.63) 


Since we are taking x_ 3> the spectrum tends to continuum, i.e. Nn —>■ N{m), 
t/^(n)(0) —>■ and the sum turns to an integral where the hrst factor in ( p.63| ) is 

the integration measure, i.e. = '^Am —>■ J dm from ( 2.57 ). Moreover, we 

can extend the integration to inhnity because the integral saturates rather quickly. 
Finally, we expand the fraction involving keeping the hrst term in the power series 
to obtain the potential: 


AKAr) = 


327rM3 


dm 


2k 


1 - 


x_ 


(2.64) 


The integral is easily calculated and the potential reads: 


= - A (1 - T) 

2 r [kry X- 


(2.65) 


where we identihed Gn = reasons to be seen later. The second part of the 

above potential is the hrst correction coming from the fact that X- is hnite. Obviously 
this correction vanishes when a;_ —> oo. Note that the above potential gives corrections 
to the Newton’s law only at distances comparable to the Planck length scale. 
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Multigravity: 4D and 5D gravity 


We turn now to the more interesting first summation in (|2.61|) in order to show that 
the coupling indeed has the “resonance”-like behaviour for Am —0 responsible for 
four dimensional Newtonian gravity at intermediate distances and the hve dimensional 
gravitational law for cosmological distances. This summation includes the KK states 
from the graviton zero mode up to the A^ 3 ;_-th level. The normalization constant in 
this region is computed by considering the series expansion of all the Bessel functions 
involved. It is calculated to be: 


JVj = 


2 4 

TT mz 


Ag{ziYx_ 




8 fc 2 

g(zi)4x_ 


( 2 . 66 ) 


where T = 4fce If we combine the above normalization constant with the unnor- 
malized wavefunction ( p.60|) , we hnd that the long distance gravitational potential is: 


Vid{r) = 




E 


nk 4:k‘^ 


327rM3 x_e^ wolziY r 

n=0 


2 I r£ , 8 fc 2 

'"'n ' 4 ' g(2i)4x_ _ 


(2.67) 


Again, since we are taking the above sum will turn to an integral with 

Z f dm. Moreover, we can safely extend the integration to inhnity since the 

integral saturates very fast for m < T/d = r~^ ke^. If we also expand the fraction 
in brackets keeping the hrst term in the power series, we hnd the potential: 


ViZ) = 


1 r , e- 

_ / (iTfl _ _ 

327rM3 io ngiziY 


2 I r 2 

m-^ + ^ 


327rM3 


dm 


32k^ 


x_7ig{ziy r [nF + 


( 2 . 68 ) 


The hrst part is the same as in the GRS model potential, whereas the second one 
is the hrst correction that comes from the fact that x_ is still hnite though very large. 
Note that the width of the “resonance” scales like something that is compatible 
with the scaling law of the masses (m„ = mrk^^), since we are working at the region 

where X- S> i.e. mn <C n'Jike~^^. The above integrals can be easily calculated in 
two interesting limits: 

For k~^ r Tc the potential is given approximately by: 


.,2x 


Vid{r < rj = 


■—( 1 -—) 

r a;_ 


(2.69) 
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Figure 2.8: The behaviour of the coupling, a{m), in the limit of X- 3> Three 
regions of interest are indicated. The region m > rus gives rise to short distance 
corrections. The mi m rUc region gives rise to 4D gravity at intermediate 
distances and 5D gravity at ultra large distances. For distances r ^ m^^, the zero 
mode gives the dominant contribution and thus we return to 4D gravity. 

where we have identihed Gn = obtain the normal four dimensional Newtonian 

potential. Note that since X- 3> the 1/x- term is indeed a small correction. 

In the other limit, r S> r^, the integrand is only significant for values of m for 
which the m^ term in the denominator of the Breit-Wigner can be dropped and the 
potential becomes hve dimensional: 


(2.70) 

7rr^ X- 

The fact that Newtonian gravity has been tested close to the present horizon size 
require that for k ~ Mpi we should have x > 45-50. Finally, we note that if we take 
the X- —>• cx) we recover the GRS result: 

\im V+__+{r,x_) = VGRs{r) (2.71) 


• Back to 4D gravity 

As we have just seen, probing larger distances than Vc, the four dimensional grav¬ 
itational potential changes to a hve dimensional one. This is the most signihcant 
characteristic of the GRS model. In the case that X- is large compared to but still 
hnite, there is another distinct region of interest, namely r ^ This follows from 
the fact that, in this limit, the spectrum is still discrete. For distances larger than of 
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the order of the corresponding wavelength of the first KK mode, the contribution to 
gravity from the KK tower is suppressed and thus the zero mode gives the dominant 
contribution, leading to the four dimensional Newtonian potential again. In this case 
the strength of the gravitational interaction is a small fraction of the strength of the 
intermediate scale four dimensional gravity. More precisely, the contribution of the 
massless graviton is l/a;_ suppressed and thus vanishes when a;_ —> oo, something that 
is expected since in this limit there is no nomalizable zero mode. The gravitational 
potential in this case is: 


V4n(r) 


1 1 
327 rM^ r 


A'oVm(o) 


r X- 


(2.72) 


Obviously this four dimensional region disappears at the limit —*■ oo since the 
spectrum becomes continuum and thus the five dimensional gravity “window” extents 
to inhnity. We should note finally that for the values of x that we consider here this 
hnal modihcation of gravity occurs at distances well above the present horizon. 


3 The Crystal Universe Model 


The last model we consider is the crystal universe introduced in [Q] (see also |l55| , [56 


It consists of an infinite array of parallel three-branes in five dimensional AdS space 
with cosmological constant A. For simplicity we assume that the branes are equidistant 
with distance, /, between two successive branes. Needless to say, in this case all the 
positive tension branes are scale equivalent. The metric ansatz that respects four 
dimensional Poincare invariance is again given by ( p.3|) where the a(y) function is 
constrained by the the Einstein equations to have the sawtooth form: 


+CXD 


(^(y) = k ^{y- ‘2jl) [29{y - 2jl) - e{y - (2j - 1)/) - e{y - {2j + 1)/)] (2.73) 


j=-oo 


The tensions of successive branes are required to be opposite and equal to ±A/k, where 
k = \j 24 m ‘^ ^ measure of the curvature of the bulk and M the five dimensional 
fundamental scale. 

We consider the general fluctuations around the vacuum ansatz of the form (|2.10| ) 
and decompose as usual the perturbation tensor h^y{x,y) as: 


h^^{x,y) = / dmh^y{m,x)'^{m,y) 


(2.74) 
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Figure 2.9: The Crystal Universe made up of an infinite array of and " branes 
with lattice spacing I and bulk curvature k. 


where T(m, ?/) is a complex function that arises from the linear combination of the 
independent graviton polarizations. This leads to the Schrodinger equation (p.l3|) for 
the wavefunction T(m, 2 ;) = \k(m, with potential: 


U(z) 


15P 


3k 

2c/(^) 


+00 


j=-oo 


Zi 


(2.75) 


The z coordinates have been dehned as usual to be the ones that make the back¬ 
ground metric conformally flat. In these coordinates the branes sit at the points 

Zj = = jzi and the function g{z) is: 


H-cxd 

g{z) = l + k {z- Z 2 j) [29{z - Z 2 j) - 0{z - Z 2 j-i) - 0{z - Z 2 j+i)] (2.76) 

j=-oo 


The solution of the Schrodinger equation for the wavefunctions for two adjacent 


cells is given [54| in terms of Hankel functions: 


4/{m,z) = 

(2.77) 

in the regions 2 ; G [ 0 , 2 ;;], 2 ; G [ 2 ;;, 22 ;;], 2 ; G [ 22 ;;, 82 ;;] and 2 ; G [ 82 ;;, 42;;] respectively. In 
the above expression is an overall normalization constant, q is the Bloch wave 



Ht ( "^giz) 
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quasi-momentum and the constants A,B,C are given by: 


A 

B 

C 


+ hfh^)ht - - 2h^hth2 

hfh^)h':^ + {hfh^ - )e2*'?^'h2 - ‘2hfh2h^ 

h~^I 12 F I 12 F 2h-^ /z2 A 
hfh2 — hi 

2htht {hth2 F hfh^)A 

^1^2 — hfh2 


(2.78) 

(2.79) 

(2.80) 


where (^^g{zi)^ and The above coefficients of the Hankel 

functions were determined by the boundary conditions, i.e. continuity of the wave- 
function at zi, 2zi, 3zi, discontinuity of its hrst derivative at zi, 3zi and the Bloch wave 
conditions relating the wavefunction at the edges of each cell^. The last remaining 
boundary condition of the discontinuity of the hrst derivative of the wavefunction at 
2zi gives us the band equation connecting q and m. We disagree with the relation 
given in 


54 and instead we hnd: 


cos{2qzi) 


{j2yi F jiy2){j2yi F jim) - yihym - 2jii2^W2 ^ 

(422/1 - 412/2) (42^1 - 4W2) ~ 


(2.81) 


where again {fg{zi)), yn = (f5'(^«)) and (f), Hn = Yn (f )• This 

dispersion relation always gives a band at zero as is to be expected intuitively. Dehning 
the parameter x = kl, we see that for 10 we may reliably approximate the width 
of the zero mode band by: 

To = 2V2ke-‘^^ (2.82) 

while the separation of the zeroth and the hrst band is 


ATi = ^ike-^ (2.83) 

with ^1 the 1st root of Ji{x). This characteristic behaviour of the above widths with 
respect to x implies that we can have a viable multigravity scenario if Tq ^ 10^®cm 
and at the same time ATj”^ ^ 0.2mm. This can happen for k ~ Mpi and x ~ 68. 

The band structure has the following form. The hrst bands have very narrow 
widths Tj with spacing ATj between them {i.e. between the i-th and the {i — l)-th 


^The Bloch wave condition reads F{zi + 2zi) = F(zi). 
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band), which are approximately: 


hi = 

Ahj = - ^i-i)ke~ 


(2.84) 

(2.85) 


where ^ 2 n+i is the (n + l)-th root of Ji(x), ^ 2 n is the n-th root of J 2 {x), and 
^ 2 n+i = f J 2 (|r+i) ’ = f jilga"*) - move on to higher bands their width is 

increasing and the spacing between them is decreasing. In the limit m —>■ cx), the 
spacing disappears, i.e. AFj —> 0, while their width tends to Fj —»• ^ke~^. In this 
limit the function f{m) in (|2.81|) tends to f{m) = cos(2z;m) and we have no forbidden 


zones. 

In order to hnd the gravitational potential we should at hrst find how the KK 
modes couple to matter and determine their spectral density. Instead of working with 
wavefunctions in an inhnite extra dimension with continuous normalization, it is more 
convenient effectively to compactify the system to a circle of length 2^/1, by assuming 
that we have a finite crystal of A/” cells [2jl, 2{j + 1)/]. We shall find the gravitational 
potential is independent of A/", so we will be able to take the limit JC —>■ oo and 
decompactify the geometry. 

In this regularization scheme we can readily see that the ever-present zero mode, 
corresponding to the four dimensional graviton, has wavefunction: 


'F(0,2:) = 


[g{z)]T^ 


( 2 . 86 ) 


and its coupling to matter on the central positive brane is: 


a(0) = 


1 


Vk 


2V2AP 2y/2M^ 


(2.87) 


The wavefunction ( p.86 ) is a smooth limit of ( 2.77 ) for m —»■ 0. Indeed, for large 


X we find the normalization constant in (|2.77|) in the region of the zero mode band is 
approximately: 

72 1 


N^, = 


zfm TaT 


( 2 . 88 ) 


The approximation of the Hankel functions with their hrst term of their power 
series proves our previous claim. Furthermore, we can see that the KK states in the 
whole zero mode band for large x will couple to matter in the central positive brane 
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as: 


a(m) 


_ 


(2.89) 


with (j){m) an unimportant m-dependent phase that does not appear in observable 
quantities. 

In the limit m ^ 00 the Hankel functions can be approximated by their asymptotic 
form. The wavefunction is then almost constant in z-space and the KK states couple 
to matter on the central positive brane as: 


a{m) 


_g* 0 M 

v(Ar4e-/2M, 


(2.90) 


where again (j){m) is an unimportant m-dependent phase. The couplings in the low 
lying bands interpolate between the above limiting cases. 

The standard procedure to calculate the spectral density is at hrst to hnd the 
number of states with masses smaller than m and then differentiate with respect to 
m. In our regularization prescription it is straightforward to count the states in the 
quasi-momentum space with masses lighter that m. The latter is simply: 


Ntotim) 


Mzi 
- Q 

TT 


SO the spectral density is its derivative: 


p(m) 


dNtotjm) 

dm 


MZi dq 
71 dm 


(T) 


Af 1 df{m) 
2 ti - f{my dm 


(2.91) 


(2.92) 


where the plus and minus signs will succeed each other for neighbouring bands, start¬ 
ing with minus for the zero mode band. For the zero band we can have a reliable 
approximation for large x of the above function: 


p{m) = Af 


vrTr 


\/l - {m/Toy 


(2.93) 


Thus, the spectral density diverges as m —Tq but the divergence is integrable 
and does not cause any problem to the following calculations. Actually, the spectral 

density diverges at the edges of every band since at these points — f{my —> 0 while 
7 ^ 0. The only point that this does not happen is at m = 0 where > 0 in 

such a way that the result is hnite. 

In the limit m ^ 00 the bands disappear and the spectrum has constant spectral 
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Figure 2.10: Spectral density for the first three bands in the crystalline model. 


density: 

p{m) = (2.94) 

kiT 

The gravitational potential, taking into account the Yukawa suppressions of the 
KK states, is simply: 

Vir) = - / dm ' \ pirn)- - (2.95) 

Jo 47r r 

For distances r 3> AF^^ we have the effective potential: 


V{r) 


1 I p{m) _ 11 

327rM2 r Jo ^ Af ^ 327i‘^M^ r Jq ^1 - ^2 

gj^l|/o(r„r) - L„(r„r)] (2.96) 


where Jq is the 0-th modified Bessel function and Lq is the 0-th modified Struve 
function. The gravitational potential is A/'-independent as expected. For distances 
AF^^ <C r <C Fq^ the above functions tend to /o(For) —> 1 and Lo(For) —> 0, so we 
recover the four dimensional Newton law with Gn = 

V(r « r„-‘) = - A (2.97) 


On the other hand, for distances r S> Fg ^ the asymptotic expansion of the difference 
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of the modified Bessel and Struve functions gives a five dimensional Newton’s law: 


V(r » r„-‘) 


2Gjv 1 

TThn 


(2.98) 


In case we take the crystal to be hnite, there will appear a region for r S> where 
gravity will turn again to four dimensional. However, as explained in the previous 
section, this region is well above the universe horizon and thus of no phenomenological 
interest. 

Finally, for distances r ^ AF^^ we will start to feel the short distance Yukawa type 
modihcations to gravity due to the presence of the bands above the zeroth one. As 

in the zero hierarchy bigravity and " H- h" multigravity model we expect these 

corrections to be important at scales of the order of the Planck length. In summary, the 
physics of the crystalline model is very similar to the GRS model as far as modihcations 
of gravity due to tensor modes is concerned. In both cases we have a crossover from 
four dimensional to hve dimensional gravity at cosmological scales. 


4 Propagator structure 


We have shown that in models in which negative tension branes are sandwiched be¬ 
tween positive tension branes there are anomalously light KK states that can lead to 
modihcations of gravity. So far we have implicitly assumed that the gravitational force 
generated by a massive graviton is identical to that generated by a massless one but 
this is not, in general, the case because the massive graviton has additional degrees 
of freedom which do not decouple in the massless limit. As a result the interaction 
generated by a massive graviton violates the normal four dimensional relation between 
the gravitational interactions of matter and light in Einsteinian gravity p^,|30[| . This 


may be seen explicitly from the form of the massive and massless graviton propagators. 

In more detail, the form of the massless graviton propagator in hat spacetime (in 
momentum space) has the form; 




2 ^ 




(2.99) 


where we have omitted terms that do not contribute when contracted with a conserved 
energy momentum tensor 

On the other hand, the four-dimensional massive graviton propagator (in momen- 
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turn space) has the form: 


_ _ 

“ 2 


1 _j_ ^ 


pi _l_ ^i 


+ ■ ■ ■ 


( 2 . 100 ) 


The difference between the two propagators 


^Qlj,v,a0 _ 


I gt^<^ga0 
6 


+ ■■■ 


( 2 . 101 ) 


is due to the additional helicity components needed for a massive graviton. There are 
three additional components needed. Two, corresponding to the graviphotons, decou¬ 
ple at low energies as they are derivatively coupled to the conserved energy-momentum 
tensor. The third, corresponding to a scalar component, does not decouple and is re¬ 
sponsible for 6G. The contribution of 6G to the one-graviton exchange amplitude for 
any two four dimensional sources and vanishes for light for which T = = 0 

but not for matter because T = Tg 7 ^ 0, hence the normal four dimensional gravity 
relation between the interactions of matter and light is violated as was hrst pointed 
out by . For example, if gravity is due to the exchange of a massive spin 2 particle 
and one insists in the validity of Newton’s law for static sources, then the deflection 
angle of light would be 25% smaller than if it corresponds to the exchange of the 
massless graviton. The fact that the bending of the light by the sun agrees with the 
prediction of Einstein’s theory to 0.02% accuracy, rules out the possibility that gravity 
is due to massive graviton exchange irrespective of how small its mass is. Of course 
there remains the possibility that a small fraction of the gravitational interactions are 
associated with a massive graviton component in the presence of a dominant massless 
graviton component. This can be realized by having an ultralight spin-2 particle with 
a very small coupling compared to graviton’s. Such an example could be the massless 
- massive graviton system in the asymmetric bigravity model case. 

However, in our discussion so far we have omitted an important gravitational mode, 
the moduli (radions) corresponding to moving negative tension branes. These radions 
as we shall see in the following chapter cancel the troublesome additional polarization 
states of the massive gravitons [p7i^,[3^,|33[| . This restores the phenomenologically 


favoured propagator structure of the above models at intermediate distances. How¬ 
ever, these helds must have unusual properties because it is known that a normal scalar 
contributes with the same sign as 6G and therefore cannot cancel it. An explicit com¬ 
putation presented in the following chapter reveals that it has a negative kinetic energy 
term leading to a contribution which cancels the troublesome 6G contribution. As a 
result the theory indeed gives rise to normal four dimensional gravity at intermediate 
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distances, but at large distances the ghost radions dominate giving in all cases scalar 
antigravity . The presence of the ghost radion held is worrying for the consistency 
of the theory. We will elaborate on this issue in the following chapter. 



Chapter 3 

Radions and dilaton 


In this chapter we are interested in describing the dynamics of the scalar gravitational 
perturbations for a system with branes in a warped background in hve dimensions. 
These excitations describe the effect of the fluctuation of the size of the extra dimension 
and/or of the relative positions of the branes. We will distinguish these two kinds of 

modes by calling the former dilaton and the latter radions [04| ,^ |^. An orbifold 
symmetry Z 2 acting on the extra dimensional coordinate as ?/ —>■ —y is also imposed. 
When the topology of the extra dimension is S^, the compact case, the Z 2 action 
has two hxed points y = 0, L. If branes are sitting on the orbifold hxed points, the 
corresponding radions are frozen. For example in a system of three three-branes where 
two of them are sitting on hxed points aX y = 0, L, we are left with just one radion 
held corresponding to the huctuation of the position of the third brane and the dilaton 
corresponding to the huctuation of the size of the orbifold. 

For the sake of simplicity we will only consider a three hat three-branes system 
embedded in AdS^ with in general diherent cosmological constants between the branes. 
This will be sufficient to clarify the role of the radion and dilaton in the multibrane 
models. We will calculate the full gravitational ehective action for this general model, 
compute the kinetic energy of the moduli and study the nature of gravity on the 
branes. As we have already noted the radion excitations play an important role in the 
context of multigravity models. They restore the correct propagator structure for the 
massive graviton excitations. However, the radions corresponding to negative tension 
branes are physical ghosts and therefore disastrous for the theory. 

Afterwards, we will concentrate on the dilaton of a hve dimensional model with 
identical four geometries and calculate its mass for bulk dynamics involving arbitrary 
minimally coupled scalar helds. A useful tool for this calculation is the sum rules 


derived by Gibbons, Kallosh and Linde ||^ which we briehy review. We hnd a generic 
formula relating the dilaton mass with the sum of the tensions of the branes and the 
curvature of the four-geometry. Demanding that this mass is not tachyonic we can 
derive the necessary condition for stabilizing the overall size of the system. 
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A2 

Ai 

Ai 

A2 






y = -L y = -r y = 0 y = r y = L 

Figure 3.1: General three three-brane model on an orbifold with different cosmological 
constants between the branes. 

1 The general three three-brane system 

We will consider a three three-brane model on an /Z 2 orbifold. Two of the branes 
sit on the orbifold fixed points y = yo = 0, y = y 2 = L respectively. A third brane is 
sandwiched in between at position y = yi = r as in Fig. p.l| . In each region between 
the branes the space is AdS^ and in general the various AdS^ regions have different 
cosmological constants Ai, A 2 . The action describing the above system is: 



where M is the five dimensional Planck mass, Vi the tensions of the gravitating branes 

and the induced metric on the branes. The orbifold symmetry y —>■ —y is imposed. 

We seek a background static solution of Einstein equations for the following four 
dimensional Poincare invariant metric ansatz: 

ds^ = a?{y)ri^j,dx^dx'^ + dy‘^ (3.2) 

The solution for the warp factor has the usual exponential form: 

{ g-fciy ,0 < y < r 

(3.3) 

g-fe2y+r(fc2-fci) ^r<y<L 


where ki and /c 2 are the curvatures of the bulk in the two regions and are related to 
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the bulk cosmological constants as: 


kl = 


Ai 


24M3 


kl 


Ag 


24M3 


(3.4) 


Moreover, the Einstein equations impose the following fine tuning between the 
brane tensions and the bulk cosmological constants: 


Eo = 24M3 ki, Vi = 2AM'' 


3(^2 - ki ) 


14 = -24M^ ko 


(3.5) 


It is straightforward to recover some special models that have been considered in 
the literature. The RSI model is obtained for ki = ^2 where the intermediate brane 

is absent (zero tension). For /ci > 0 and /c 2 = —ki we get the " H- V” bigravity 

model |2l|,|5^. For /ci > 0 and /c 2 > ki we get the " + H —" brane model considered 
In the decompactification limit where L ^ 00 we get also two interesting 


m 


26 


models: For ki > 0 and /c 2 = 0 we obtain the GRS model |^| and for fci > 0 and 
/c 2 = 0 > /ci the non-zero tension version [R^ of the model considered in p7 . 


2 Effective action 

Our purpose is to study fluctuations of the background (|3.2|). The first important 
observation is that there exists a generalization of Gaussian normal coordinates such 
that in the perturbed geometry the g^y components of the metric are zero and at the 
same time the embedding of branes is still described hj y = 0, y = r and y = L (see 
for instance the Appendix of [^). In that case the metric takes the form: 

ds‘^ = g^y{x,y)dx^dx’' + gyy{x,y)dy‘^ . (3.6) 

When analyzed from a four dimensional point of view, in each region, perturbations 
are of of three types. 

• Spin two: 

Tensor-like perturbation hy^^{x,y) corresponding to massive (and massless) four 
dimensional graviton^: 




a^iy) 


+ h^u{,x,y) 


dx^dx'^ -I- dy^ 


(3.7) 


^We denote hy^i, the conformal perturbation in contrast to the previous chapter graviton pertur¬ 
bation hy^. Their relation is hyi, = a^hy^. 
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• Spin zero: Dilaton 


Scalar perturbation fi{x) corresponding to the fluctuation of the size of the 
orbifold 


^ Q{y)fi{x)] 7]^^ dx^dx'" + [1 + q{y)fi{x)] dy^ (3.8) 


• Spin zero: Radion 


Scalar perturbation f 2 {x) corresponding to a fluctuating distance of the inter¬ 
mediate brane 


ds'^ = a^(y) [(1 + B(y)f 2 {x)) + 2e{;y) d^d„f 2 {x)] dx'^dxd' + [1 + A{y) / 2 (x)] 

(3.9) 


In the above decomposition we should note that we have kept only the zero modes 
of the radion and the dilaton. Their massive KK-states are “eaten” by the graviton’s 
corresponding modes in order to give them mass, in a geometrical realization of the 
Higgs phenomenon. The function e{y) is needed in order that the metric satisfies the 
Israel junction conditions in the presence of moving branes. As it will be shown later, 
the values of dye at ?/ = 0, r, L are gauge invariant. When ki ^ k 2 a non-trivial e is 
required. 

The generic perturbation can be written as: 


=a?{y) |[1 ^i{x,y)] + 2e{y) dyd^f 2 {x) + hyu{x,y)^ dx^dx'" 

+ [I + ^p2{x,y)]dy‘^ 


(3.10) 


where 


Vi{.x,y) = Q{y)fi{x) + B{y)f 2 {x) 
^ 2 {x,y) = q{y)fi{x) + A{y)f 2 {x) 


(3.11) 


Given the expression (|3.3|) for a, the Israel junctions conditions (continuity of metric 
components and jump conditions) at y = 0,r, L, simply require that A, B, dye, Q, q 
are continuous there [^. The four dimensional effective action Sefj for the various 
modes is obtained inserting the ansatz (|3.11D in the action (|3.1| ) and integrating out 
y. So far the functions A, B, e, Q, q have not been specified, however imposing that 
Sgff contains no mixing terms among h and fi one determines Q, q and finds A in 
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terms of B (see Appendix A). The latter satisfies: 


A. 

dy 





(3.12) 

/ da\ dB 
\dy) dy ° 

(3.13) 


The condition that h and /j do not mix between them allows us to give a clear 
particle interpretation for the quanta of h and /j. However, despite the fact that we 
have chosen that gauge, the results that we obtain are bound to be the same as in the 
general case where there would be quadratic mixing between the modes, because of 
the five dimensional gauge invariance of the theory. As a consequence of the no-mixing 
conditions the linearized Einstein equations for (|3.11|) will consist in a set independent 
equations for the graviton and the scalars. 

The effective Lagrangian reads (see Appendix A): 


Jd‘xC,„ J d^X {Cora. + 

Ca„ = 2M=> j dy (A) + j sjr 


(3.14) 


^Scal 


with 


^Scal /Ci/in/i + /C2/2n/2 

/Cl = 2M3 -c2 r a -2 dy, IC^ = -2M^ ^ ± dy 

2 J-L iJ-L \dyj dy' 


(3.15) 


where c is an integration constant. 

In (|3.14|) , the spin-2 part, as expected, contains the four dimensional Pauli-Fierz 

Lagrangian cAfA) for the gravitons without the mass term (see Appendix A), plus a 
mass term coming from the dimensional reduction. In the scalar part Cscai the mass 
terms are zero since fi are moduli fields. The indices in the above expressions are 
raised and lowered with 77 ^,^. Notice that after the no-mix conditions are enforced, the 
effective Lagrangian contains the undetermined function e{y). 

The metric ansatz in (p.ll|) is related to a special coordinate choice, never¬ 
theless, a residual gauge (coordinate) invariance is still present. Consider the class 
of inhnitesimal coordinate transformations —*• x'^ = x^ + such that the 

transformed metric = Gj^j^ + 5G retains the original form ( p.lOj ) up to a redef¬ 
inition of the functions q,Q, A, B,e and the dilaton and the radion held. Consistency 
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with the orbifold geometry and the requirement that the brane in ?/ = r is kept fixed 
by the diffeomorphism lead to 0) = r) = L) = 0. From 

■sGm* = (a.w) 

one can show that has to be of the form: 


^^{x,y) = i^{x) -W{y)r]f^''d^f 2 {x), ^^{x,y) = a^W'{y) f 2 {x) 

with hF'(O) = W'{r) = W'{L) = 0 


(3.17) 


The case VF = 0 corresponds the familiar pure four dimensional diffeomorphisms, 
under which transforms as spin two field, /, as scalars and q,Q, A, B,e are left 

unchanged. On the contrary the case = 0, W 7 ^ 0 is relic of five dimensional 
diffeomorphisms and one can check that q, Q, A, B, e are not invariant and in particular 
6e = W. As a result, the values of dye in 0, r and L are gauge independent and this 
renders Ceff free from any gauge ambiguities. 


3 Scalars Kinetic Energy 


3.1 The compact case 

In this section we will focus on the part of the effective Lagrangian involving the 
scalars and concentrate on the dilaton and radion kinetic term coefficients /Ci and 1 C 2 - 
In particular we are interested in the cases when the radion becomes a ghost field, i. e. 
IC 2 < 0. Firstly, it is trivial to obtain the dilaton kinetic term /Ci by carrying out the 
first integral of (|3.15|): 


/Cl = 


a ^(r) — 1 a ‘^{L) — a ^(r) 
h /C 2 


(3.18) 


It turns out that for any possible values of fci, ^2 and r, L the above quantity 
is positive definite. The radion kinetic term on the other hand is more involved. 
Integrating (|3.12D we get the radion wavefunction for the two regions (y > 0): 


{ Cl a ^ + 2ki dye , 0 < y < r 

(3,19) 

C 2 a ^ + 2k2 dye ,r < y < L 


where ci and C 2 are integration constants. The orbifold boundary conditions demand 
that dye{0) = dye{L) = 0 since e is an even function of y. From the non-mixing con- 
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ditions for radion and dilaton ( |3.13|) and the continuity of B we are able to determine 
C 2 and dycir) as the following: 


^2 a?{r) — 1 

C2 = Cl — 


ki ( a(r) \ . 

l^a(L) J 


e r = 


Ci/C2 


2 / ci (/ c 2 - ki)a‘^{r) 


ki a^{r) — 1 

k2 


a(L) 


- 1 


(3.20) 


(3.21) 


Therefore, the values of the radion wavefunction B at the brane positions are given 
by the following expressions: 


B{r) = 


5(0) = Cl 

Ci/C2 1 — 0?{,L) 


{k2 - ki) 


a2(L) 


( £(!i^ _ 1 


BiL) = 


Ci/c 2 a^(r) — 1 


ki 


a?{L) 


ah) 

a(L) 


(3.22) 

(3.23) 

(3.24) 


Thus, we can carry out the second integral on ( p.l5|) to hnd the radion kinetic term 
coefficient: 

/C 2 = 3M3 1 - 1") B\r) a\r) + ^ B\L) a\L) - ^ B^{0) 

/C2 / fv2 


3M^cl 

ki 


k 2 a2(r)(a^(L) - 1)2 k 2 {a^{r)-l) 


{k2 - ki) 


a\L) 


H 2 


ah) 

a(L) 


- 1 


+ 


ki 


a\L) 


a(i-) 

a(L) 


- 1 


(3.25) 


The above quantity is not positive dehnite. In particular, it turns out that it is 
positive whenever the intermediate brane has positive tension and negative whenever 
the intermediate brane has negative tension. This result is graphically represented in 
Fig.^]^ where the (fci, ^ 2 ) plane is divided in two regions. 


3.2 The non-compact limit 

It is instructive to discuss the decompactification limit in which the third brane is sent 
to inhnity, i.e. L —> +cx). To examine this limit we distinguish two cases: 
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/Co >0 


C 


B' / 
/ 

^ A' 


^ ki = k2 


A 


B 


C GRS 


/Co < 0 


ki 


Figure 3.2: Sign of the radion kinetic term in the {ki, /C 2 ) plane. In the regions A, B', 
C the moving brane has positive tension and the radion positive kinetic energy. In 
the regions B, C, A the moving brane has negative tension and the radion negative 
kinetic energy. We show the GRS line for the non-compact case. The dashed line 
corresponds to ki = ^ 2 , i.e. a tensionless moving brane. 


The case ko > 0 


In this case we have a(oo) cx: lim e 

L^oo 


-k2L _ 


= 0. The dilaton kinetic term is trivial since 


from (|3.18|) we obtain /Ci —> cx). In other words the dilaton decouples from the four 
dimensional effective theory and the condition of absence of kinetic mixing between 
the scalars (|3.13| ) plays no role. The radion kinetic term coefficient can be read off 
from (|3.25|) : 


/C2 = 


3M^cl 

ki 


^2kir 


ko 


(/C2 - ki) 


- 1 


(3.26) 


This result agrees with the computation of |^. It is easy to see that the radion 
has positive kinetic energy when the moving brane has positive tension and negative 
kinetic energy when the tension is negative. Indeed, for 0 < /ci < /c 2 , or for /c 2 > 0 
and /ci < 0 we have a positive brane and positive kinetic energy. On the other hand, 
for 0 < k 2 < ki we have a negative tension brane and negative kinetic energy. In the 
limit /c 2 —*■ 0 we get the GRS model with negative kinetic energy as in |32 . 


The case k 2 < 0 


In this case a{oo) oc lim e > 00 . This time the dilaton plays in the game since 

L—^oo 
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/Cl is finite as seen from (|3.18|) and has the value: 


/Ci = 


3M3 c2 

2/Ci|/C2| 


[e^'^^^{\k2\ + h)-\h\\ 


(3.27) 


which is manifestly positive dehnite. The presence of dilaton mode is somewhat sur¬ 
prising since it describes the fluctuations of the overall size of the system which is 
inhnite. Something similar happens in the Karch-Randall (KR) model that has 
a dilaton mode although it is non-compact. The dilaton is a remnant of the decom- 
pactihcation process of the " -|- +" model which we will examine in the following 
chapter and enters in the game because the inverse of the warp factor is normalizable. 

The radion kinetic term coefficient can be obtained from taking the limit L —> -|-cx) 


in (|3.25|). We get 


/C2 = 


3M3, 

ki 


—2k\r 


ko 


(/C2 - ki) 


(3.28) 


The same considerations for the compact case applies here. When fci < ^2 < 0 we 
have a positive tension brane and a positive kinetic energy. On the other hand when 
k 2 < ki < 0, or for /c 2 < 0 and /ci > 0 we have a negative tension brane and a negative 
kinetic energy. In the GRS limit /c 2 —^ 0, the radion has negative kinetic energy. 


4 Gravity on the branes 

In this section we will study how the moduli and graviton(s) couple to matter conhned 
on the branes. For simplicity we will study the non-compact case. For this purpose, 
we consider a matter Lagrangian where we denote generically with <Fj 

matter helds living on the branes; is the induced metric on the i-th brane. Our 
starting point is the following action: 

S= J d^xdy [2M^R - A( 2 /)] j d^xsj-G^Wi 

* (3.29) 

+ Y, j d^xdy\l^)Cm{^j,G^^) 

i 

We have already calculated the effective action for the gravity sector for the per¬ 
turbation (^.101) in the previous chapter. It is useful to decompose the perturbation 
h{x,z) in terms of a complete set of eigenfunctions^ of the graviton kinetic 

^We denote by 'I'(z) the gravitons wavefunctions in the conformal frame in contrast to the previous 
chapter graviton wavefunction 'I'(z). Their relation is 'I'f^) = a^4>{z). 
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operator (a suitable four dimensional gauge fixing like de Bonder is understood): 


Hx,z) = + / dm^{y,m)h^^\x) 


(3.30) 


having taken into account both the discrete and the continuum part of the spectrum. 
The effective four dimensional Lagrangian has the following form: 


C = + 2M3 ^ CpF{h^^\x)) + /Cl /iD/i + /C2 h^h 






S(9.) 


(3.31) 


The matter helds have been rescaled to make them canonically normal¬ 

ized and the energy momentum tensor for matter localized on the Tth brane is 
dehned with respect to the rescaled helds <h^. By construction the induced background 

metric on the branes is the hat four dimensional Minkowski metric. The asterix de¬ 
notes that the sum has to be converted into an integral for the continuum part of the 
spectrum. Finally, dehning the canonically normalized helds: 

= /, = y^/i. /2 = y «/2 (3.32) 


the Lagrangian reads: 


£ = £„(<!.;,.,) + £pf(S<">(i)) + i /,□/, + i sgn(/C2) huh 




2" " 2 

Q{yt) 




(z) StlTTli 


B{y^) 

2v^ 


(3.33) 




Thus, the dilaton, the radion and the graviton (whenever it is normalizable) have 
the following couplings respectively: 


Cd — 


Q{y. 


br ) 




Cr — 


B{ybr) 


Cg — 




br ) 


2V2M^ 


(3.34) 


For convenience we will dehne Cn = and write the above quantities 


as: 


(b 

Mi) _ p 9 d 

, 


(b 

p{i) _ p 9 r 


Cg — Cn- 


¥^\ybr) 

VT^ 


(3.35) 
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Figure 3.3: The function a{y) = — log[a(|/)] for all possible combinations of ki, /c 2 - 
The regions are named in accordance with the Fig.^]^ phase diagram. 


where the dimensionless radion couplings are: 



9i-k2) 


2kirtl + \hA — 1 
^ \k2\ ^ 



1 


p2fci rsgn(fc9~) fc 2 _ 1 

^ k2-ki 



2fcir (0) 

Ud 


(3.36) 




k2 

k2 - ki 



(3.37) 


We should note at this point that there are certain cases as we will see in the 
following that multigravity is realized and four dimensional gravity in intermediate 
distances is due to more than one mode. In that case we will denote by Cq the 
coupling of the “effective zero mode” even though there might not be a genuine zero 
mode at all. 

We will now discuss how these couplings behave is the six distinct combinations of 
ki and /c 2 shown in Fig.|3.3|. 


Region A 


In this case, we have the non-compact system of two positive tension branes dis¬ 
cussed in [^,^,2^. The volume of the extra dimension is hnite and therefore we have 
a normalizable graviton zero mode. The KK tower will be continuum and its coupling 
to matter on either branes will be exponentially suppressed for reasonably large r as 
can be easily seen by the analysis of |2^. The graviton wavefunction in the conformal 
gauge is constant, i.e. T^°^(|/) = A. Thus the universal coupling of the zero mode will 
be Cg = Cn- 
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The dilaton in this region will simply be absent, or equivalently it will have zero 
coupling Cd = 0 to the branes because JCi —*■ oo. The radion in this region has positive 
kinetic term and from (|3.37|) we find that its coupling to matter on the y = 0 and 
y = r branes is bounded as following (see Fig. |T^) : 




(0) ^ -fcir 


Sr ^ e*" 


(3.38) 


The radion on the central brane is always weakly coupled and decouples in the 
tensionless moving brane limit. On the other hand the radion is always strongly 
coupled on the moving brane and diverges in the tensionless limit where ^2 —*■ /ci. 


Region B 

The volume of the extra dimension is still finite and therefore we have a normal¬ 
izable graviton zero mode. When ki = k 2 the brane is tensionless and as soon as 

^2 < ki the KK continuum starts to develop a “resonance”. This “resonance” is 
initially rather broad as /c 2 —> ki and its width decreases until it coincides with the 
GRS width T ~ /ci exp(—3/cir) as /c 2 —> 0. In this region four dimensional gravity at 
intermediate distances is the net effect of the massless graviton and the lower part of 
the KK continuum that contributes more and more as ^2 0. 

The dilaton is again absent, Cd = 0, since still we have /Ci —> cx). The radion in 
this region is a ghost field, to compensate for the presence of the extra polarization 
states of the contributing massive gravitons, and its coupling to the branes is bounded 
as (see Fig.U)- 

0 < < 1, < 0 (3.39) 

The radion on the central brane interpolates between the decoupling limit of the 
tensionless moving brane to the GRS limit where it couples with strength equal to 
the one of the graviton. On the other hand, the radion coupling on the moving 
brane interpolates between the infinitely strong limit of the tensionless brane (strictly 
speaking there is no brane in y = r and we are left with RS2) to the GRS limit where 
it decouples. 


Region C 

In this region we have a infinite analogue of the "H-model with no normalizable 

zero mode present. The graviton spectrum for /c 2 —^ 0 is approximately equally spaced 
and there is a resonance which coincides with the GRS width and gets more and more 
broad as /c 2 gets more negative. Soon enough a special light state is singled out as 
in the " -I - \-" model, a behaviour that persists for all values of k 2 in this region. 
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From the above behaviour we deduce that four dimensional gravity at intermediate 
distances is generated by the lower part of the discrete spectrum as /c 2 —*■ 0, whereas 
only by the special state for all other values of ^ 2 - 

The dilaton in this region is present and its coupling to the branes is bounded as 
following (see Fig.^.4|): 


0 < (3.40) 

It is always weakly coupled to matter on the central brane and becomes strongly 
coupled but saturated on the moving brane. 

The radion in this region is again a ghost held to cancel the unwanted extra massive 
graviton polarization states and its coupling to the branes is bounded as (see 

gf~l, 0 < gj' < e*"’- (3.41) 

On the central brane it couples always with strength equal to the one of the “effec¬ 
tive graviton” and on the moving brane it interpolates between the decoupling limit 
of the GRS case to a strongly coupled region with saturated coupling as the tension 
of the second brane gets inhnite. 


Region A' 


In this region we have no normalizable zero mode and the two negative tension 
brane system resembles an inverted version of the RS2 model. The KK spectrum is 
discrete and all excitations he above the characteristic curvature scale ki. Thus the 
low energy effective theory does not have four dimensional gravity at all. 

The dilaton in this system is present and the coupling to the branes is approxi¬ 


mately constant (see Fig. |3^) : 


~ 1, gV-I (3.42) 

The radion held is a ghost and its coupling to matter on the branes is bounded as 
(see Fig.|3^): 

gf > , 0 < gj> < (3.43) 

The radion on the central negative tension brane is always weakly coupled and 
vanishes in the limit of /c 2 —*■ /ci = — |/ci | of the tensionless moving brane. On the other 
hand the coupling on the moving brane is bounded from below in the limit of inhnite 
negative tension brane and diverges at the limit of the tensionless brane. 
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Figure 3.4: The dimensionless couplings of the radion (upper) and the dilaton (lower) 
to matter on the y = 0 or the y = r branes. The left diagrams correspond to fci > 0 
while the right ones for ki < 0. The regions A, B,... are in accordance with the ones 
of the phase diagram in Fig.p.2|. The diagrams are not in scale. 


Region B' 

In this region there is again no normalizable zero mode and the system of the 
negative and positive tension branes still resembles the inverted RS2 model. The 
KK spectrum is almost identical with the one of the previous region except for the 
limit k 2 —>■ 0 when the spectrum drops below the curvature scale ki and a continuum 
develops. 

The dilaton coupling to the branes decreases from the constant value of the previous 
region, to zero as /c 2 —»■ 0 (see Fig.|3.4|): 

0<^7g^<l, 0 < (3.44) 


The radion has positive kinetic energy and its coupling to the branes is bounded 
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as following (see Fig.|3.4|): 


0 < < 1, < 0 (3.45) 

It is weakly coupled in the central brane and the coupling interpolates between 
zero from the previous region to one in the inverted GRS case where ^2 —^ 0. On the 
second brane it is divergent as ^2 —> fci = — |A;i| and vanishes as /c 2 —^ 0. 

Region C 

In this region we again have a system of a negative and a positive tension brane, 
but gravity can be localized on the moving positive tension brane. There is a normal¬ 
izable zero mode that mediates four dimensional gravity and a KK continuum with 
suppressed couplings on the branes. 

The dilaton in this region will be absent, or equivalently it will have zero coupling 
Cd = 0 to the branes because K-i —>• cx). 

The radion will have again positive kinetic term and its coupling will be bounded 
as following (see Fig. |3.4|) : 


(3.46) 

It is approximately constant on the central brane with strength equal to the one of 
the graviton and is weakly coupled to the moving brane with coupling interpolating 
between zero in the inverted GRS case with /c 2 0 to a certain value as the tension 

of the moving brane becomes inhnitely large. 


5 Review of sum rules 


In this section we will review the sum rules presented in in order to discuss the 
dilaton stabilization in the presence of bulk dynamics of arbitrary minimally coupled 
bulk scalar helds. We will concentrate as in in the case where the background 
metric of the hve dimensional spacetime can be written in the form; 


ds'^ = a{yYg^j,^{x)dx^dx'' -\- dy'^ 


(3.47) 


with g^v{x) a general background four dimensional metric and a{y) a generic warp 
factor. We should stress here that this is not the most general choice of metric in hve 
dimensions as we have explicitly assumed that the all four dimensional sections have 
the same conformal geometry. The particular choice g^y = is valid for the models 
which we have discussed so far. 
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We can now consider an arbitrary number of minimally coupled scalar bulk fields 
with internal metric Qjj and arbitrary bulk potential ld($) (which includes 
bulk cosmological constant), coupled to an again arbitrary number of branes with brane 
potential Aj(<f)) (which again includes the brane tensions). The action describing the 
above system is the following: 


S = J d^xdy^ (^2M^R - - f"($) - ^ A4<h)5(|/ - 

(3.48) 

where is the induced metric on the brane and M the fundamental hve dimensional 
scale. The Einstein equations arising from the above metric can be written in the form: 


= 

1 4 

rppL rp^ 

3 35 

(3.49) 

4M^Rl = 

-Ir; + |r| 

(3.50) 


where the energy-momentum tensor components are: 


T; = -a^<h-a'^<h-24>'-$'-4E(<h)-4 5^Ai($)5(|/-|/i) (3.51) 

i 

r| = ■ $'-1/($) (3.52) 

with the indices in the above formulas raised and lowered by = a{yYg^y{x) and 

where dot product denotes contraction with the internal metric Qij. Since we are 

interested in a background conhguration, the <9^$ ■ terms can be dropped. The 
Ricci tensor is easily calculated to be: 

R/) = a-‘^Rg-l2{a'fa-‘^-Aa''a-^ (3.53) 

Rl = -4a"a-^ (3.54) 


If we now consider the function (a”)" with n an arbitrary real number, its integral 
around the compact extra dimension is zero. Using ( p.49| ), (|3.50|) , ( p.53|) , (|3.54|) we 
arrive at an inhnite number of constraints [4T]|: 


dy + {2n - 4)T|) = - n)Rg (p dy a 


^n—2 


(3.55) 


As it is obvious, these constraints are a natural consequence of the equations of 
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motion and are totally equivalent to them. It is straightforward to see for example 
that they are satished in the RSI model |]^ and the bigravity/multigravity models 
pT|,p6|,^. We will single out three constraints which we will be important for the 
subsequent discussion, namely the ones for u = 0,1, 2: 


dy (TH - AT^) = AM^Rg (b dy a 


dya{T!b-2T^) = 0 


dy 


-2 


-AM^Rg (f) dy 


(3.56) 

(3.57) 

(3.58) 


One could also use the above constraints for non-compact models, but should be 
careful that the above derivation makes sense. For the RS2 model , all constraints 
are valid for n > 0. For the GRS model only the n = 0 constraint is valid and for 
the KR model [Q (see next chapter) all constraints are valid for n < 0. 


6 Dilaton in warped backgrounds 


We now consider the perturbation related to the overall size of the compact sys¬ 
tem, namely the dilaton. As previously noted in the chapter, we do not consider the 
KK-tower of the dilaton since these states can be gauged away by hve dimensional 
diffeomorphisms. The general form of the metric for the physical radion perturbations 
that do not mix with the graviton(s) is given in P2| , [33[1 . For the dilaton, the ansatz is 
rather simple and following [Q we write it in the form^ 

ds^ = ^'^^^^a{yYggy{x)dx^dx^ + (^ + dy^ (3.59) 


Substituting the above metric in the action (|3.48|) (see Appendix B for analytic 
formulas), integrating out total derivatives, throwing out y-independent parts and 
keeping terms up to quadratic order, we get: 


S= d‘^xdy^/^ 


^We are interested only in the quadratic in ^{x) Lagrangian derived from the above metric, so 
only the part of the metric up to linear order in ^( x) is i mportant. However, we prefer to keep the 
same notation as and do the calculation with (|3.59| ). The result, of course, for the quadratic 
Lagrangian would be the same if we used the metric up to linear order in 7 ( 0 ;). 




(3.60) 
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with 


£1 = 2M3(4(a')^ + 16a"a) + ^a^<l>'-$' + aV(<h) + 2a2^Ai(<l>)(5(2/-2/0 (3.61) 

i 

C 2 = 2M^{a-‘^Rg-32{a'fa-^ + 32a''a-^) + W-^' + A^\i{^)5{y-yi) (3.62) 


At this point, let us work out the integral over the extra dimension of the tadpole 
term £1 of the Lagrangian. This gives: 


dyCi = —- (p dy a 


- 2T| - 16M=^ ( ^ + 4— 


(3.63) 


where we used the energy-momentum tensor components found in the previous section 
with respect to the unperturbed background metric (|3.51|) , (|3.52|) . We can further 
simplify this quantity if we use the equations (p.49|) , (p.50|) , (|3.53|) , (|3.54|) which hold 
for the background metric. The resulting expression is: 


dyCi = -(fdya^ [t; + 4M^a-^Rg] 


(3.64) 


which is exactly zero because of the n = 2 constraint (p.58 ). This result should 
have been expected since the perturbation (p.59|) is bound to extremize the effective 
potential when one evaluates the action using the background equations of motion. It 
is worth mentioning here that in the case that the warp factor is effectively constant 
(a ~ 1), as it was assumed in ||^, the condition that the expression (|3.63|) vanishes, 
is identical with the n = 1 constraint ( p.57|) . 

Our next task is to read off the mass of the dilaton from the action functional. For 
this reason we define the canonically normalized dilaton field with mass dimension one 

7 ^ = (6M^ ^ dya~^) 7 ^ = ^ 7 ^. Then the mass of the canonical dilaton 7 is: 


(3.65) 


dyR2 

After a lot of simplifications using the relations ( p.4g|) -( p.54| ) we obtain: 


m? = 


^ ^ dy{10M^a-^Rg + $'■$' + 4 ^ A,($)5(|/ - y^)) (3.66) 


We can further simplify the expression using the n = 0 constraint (|3.56|) and get a 
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more suggestive result: 


= - 2M^a-^Rg) 


(3.67) 


or equivalently, 


= 




(3.68) 


From the second expression it is clear that we cannot have a massive dilaton if the 
sum of the effective tensions of the branes Aj(<F) is exactly zero and at the same time 
they are kept flat. This is the same conclusion that appeared in regarding the 
Goldberger-Wise mechanism [M in the RSI scenario. Regarding the latter mecha¬ 
nism this is not a surprise, since promoting the Goldberger-Wise to an exact solution 
means that either the scalar held energy will generate a four dimensional cosmological 
constant, or if one insists on having hat branes, this will introduce a new and diher- 
ent hne tuning of the tensions which gives a mismatch of brane tensions. Moreover, 
the absence of a tachyonic mass will guarantee the stabilization of the overall size of 
any system with the above characteristics. By (|3.67| ), ( |3.68| ) we get two equivalent 
conditions: 


d|/(<F'■ $'- 2MV2Rg) > 0 


J2Xii^)+3M^Rg ^ 


dya 


-2 


< 0 


(3.69) 


(3.70) 


If one wishes to have hat branes, then the sum of brane tensions should be negative 
or equivalently one should have a non constant (in y) scalar held conhguration. In the 
case that the above expressions (and thus the mass) vanish, one should look for the 
higher orders of the ehective potential to examine the stability of the system. 

Finally, we should note that the formulas (|3.67|), (|3.68|) are valid even for non¬ 


compact models whenever the dilaton mode is normalizable. This happens for example 


in the KR model and one can hnd from the above expressions the mass of the 
dilaton, in agreement with [^]. In that case, the dilaton mode cannot be attributed 
to the fluctuation of the overall size of the system, but can be understood to be a 
remnant mode if we start with the compact " -|- system and send one of the 
branes to inhnity (he. after decompactifying the system). We will see this possibility 
in the next chapter. 


























Chapter 4 

AdS/i^ brane multigravity models in 
five dimensions 


As we have seen clearly in the previous two chapters, the flat brane multigravity 
models in flve dimensions necessarily have moving negative tension branes. These 
in turn have corresponding radions with negative kinetic energy which are disastrous 
for the theory. It is clearly desirable to have a model with only positive tension 
branes. In this chapter we will demonstrate how one can obtain such a model in flve 
dimensions. The price that one has to pay is that the branes are no longer flat, but 
instead AdS^. A two brane model with only positive tension branes was already known 

and exhibited a characteristic bounce of the warp factor. We 


m 


know from chapter ^ that a model that has a bounce in the warp factor is bound to 
have bigravity, something that we will explicitly show for this particular model in the 
present chapter. Furthermore, we can construct a crystal model with an infinite array 
of positive tension branes, but we cannot reproduce the quasi-localized model because 
the presence of the cosmological constant on the branes prevents the warp factor from 
being asymptotically constant. 

Since there is no ghost modulus in these configurations one may worry about the 
structure of the graviton propagator at intermediate distances. As we will show ex¬ 
plicitly, the AdS 4 ^ brane configurations can evade this problem because of an unusual 
property of the graviton propagator in AdS background. We will show in the second 
part of this chapter that the van Dam-Veltman-Zakharov discontinuity can in certain 
cases be circumvented and a particular example is the the " + +” model. In more 
detail, the troublesome additional polarizations are very weakly coupled and do not 
disturb the normal structure of the graviton propagator. 

Nevertheless, the construction that we obtain faces phenomenological difficulties 
because of the remnant negative cosmological constant on the brane. This sets an 
horizon scale and unfortunately the Compton wavelength of the light state of the sys¬ 
tem lies exponentially far from this horizon. This does not change even if we consider 
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Figure 4.1: The " + +" bigravity model with two AdS^ "+" branes at the hxed points 
of a S^/Z 2 . 

a highly asymmetric version of this model. As a result, this bigravity model makes 
no predictions for observable deviations from Newtonian gravity at ultra-large dis¬ 
tances. Furthermore, the very fact that the remnant cosmological constant is negative 
is observationally disfavoured. 

1 The " + bigravity model 

The model consists of two three-branes with tensions Vi and V 2 respectively, in an 
AdS^ space with five dimensional cosmological constant A < 0. The hfth dimension 
has the geometry of an orbifold and the branes are located at its hxed points, i.e. 
Lq = 0 and Li = L (see Fig.[4.1|). Due to orbifolding, we can restrict ourselves to the 
region 0 < y < L. Firstly, we hnd a suitable vacuum solution. The action of this 
setup is: 


Sg= [ d^x [ dyV^[2M^R - A] - V [ d^xVi\f^( 

J J-L ,• Jy=Li 


(4.1) 


where C’fH is the induced metric on the branes. The Einstein equations that arise from 
this action are: 


Rmn - -GmnR - 


AG 


MN 


J-Gb) 


(4.2) 


In order to hnd a specihc form of the equations of motion we need to write a metric 
ansatz which will take in account the spacetime symmetries of the three-brane. Since 
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we would like not to restrict our model to flat solutions on the branes, we should make 
a choice which will let us interpolate between the maximally symmetric space-times 
in four dimensions, i.e. de-Sitter, Minkowski and Anti-de-Sitter. The metric ansatz 
that accomplishes this is the following: 


ds^ 


a\y) 


rj^jydx^dx’' + dy"^ 


(4.3) 


where x"^ = rj^^x^x^. It can be shown that the Ricci scalar for this metric is i? = 
yl2H‘^. From now on, the upper sign corresponds to AdS 4 space while the lower to 
dS 4 space and always > 0. Thus, this metric represents all maximally symmetric 
spaces. A straightforward calculation of the Einstein’s equations gives us the following 
differential equations for a{y): 


(a'fo))^ = (4.4) 


By solving the above equations we And that the solution can be written in the 
form: 


with 


= 


Vik 

a{y) = cosh(fc||/|) — sinh(/c||/|) 


- A2) , ^ < El for dS 4 branes 

0 , = Ri for flat branes 

fcRA2 4 ^ 2 u 2 \ |A| 


(4.6) 


(4.7) 


^(A^ — Vik"^) ,-Y>Vi for AdS 4 branes 
where we have normalized a(0) = 1 and assumed R > 0. Also we have defined 


k = 


- / -A 

24M3 • 


Additionally, in order to have this solution, the brane tensions R, R2, the bulk 
cosmological constant |A| and the position of the second brane L must be related 
through the equation: 

R + R 


tanh(/cL) = k\A\ 


|A|2 + PRR 


(4.8) 


Let us now restrict ourselves to the case of AdS 4 spacetime on the two branes which 
will turn out to be the most interesting. In this case the condition Ri must hold. 
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Hence, we can define tanh(/c|/o) = and write the solution in the form: 


a{y) = 


cosh(fc(|/o - \y\)) 

cosh(%o) 


(4.9) 


from which it is clear that the warp factor has a minimum at y = yo- From this point 
we can see the role of the AdSi on the branes, i.e. the role of the condition ^ > Hi. 
This condition allows us to have the bounce form of the warp factor and the second 
brane to have positive tension and give us, as we will see shortly, a phenomenology 

quite similar to the " H- [-" bigravity model [^]. This can be easily seen from the 

(|4.8|) which relates the brane tensions and the distance between the branes. From 
this equation we indeed see that by placing the second brane after the minimum 
of the warp factor we can make the tension of the second brane positive and thus 
both branes that appear in the model have positive tension avoiding the problems 
associated with negative tension branes. In fact it is clear that the present model 

mimics the characteristics of the " H- V" bigravity model since what we effectively do 

is to reproduce the effect of the presence of a negative tension brane, i. e. the bounce 
form of the warp factor, with another mechanism allowing a negative four dimensional 

cosmological constant on the brane. Note that in the limit that Hi —> (fiat limit) 
the minimum of the warp factor tends to infinity and if we wish to have a second brane 
at some hnite distance, it will necessarily have negative tension. 

The relationship between the four dimensional effective fundamental scale M* and 
the five dimensional fundamental scale M can be easily found by dimensional reduction 
to be: 


m! = 


k cosh^{kyo] 


[kL + sinh(/cL) cosh(fc(L — 2yo))] 


(4.10) 


The above formula tells us that for finite L the compactification volume is hnite 

and thus the zero mode is normalizable. In the case where we send the second brane 
to inhnity, the compactihcation volume becomes inhnite which indicates that the zero 

mode becomes non-normalizable. Note that M* is not necessarily equal to the four 
dimensional Planck scale Mpi since as will see shortly, at least for a sector of the 
parameter space of our model, gravity is the result not only of the massless graviton 
but also of an ultralight KK state. 


The warp factor renormalizes the physical scales of the theory as in . Thus, all 
mass parameters mo on a brane placed at the point ybr are rescaled as: 


m = a{ybr)mQ 


(4.11) 
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Hence, one can choose a distance between the two branes such that this rescaling 
leads to the creation of a desired mass hierarchy. Let us note at this point that 
this interbrane distance, and thus the resulting hierarchy, is dictated by the relative 
strength of the tensions of the two branes. This is not the case as the flat brane models 
where the hierarchy is to be hxed by some stabilization mechanism. The AdS brane 
system is self-stabilized and this can be easily seen by the fact that the dilaton in this 
case is massive. Indeed from ( p. 68 ) we obtain = 4H^ > 0 which means that the 
system is truly stabilized. 

However, since we consider non-flat solutions on the branes, we have to make sure 
that the four dimensional effective cosmological constant does not contradict present 
experimental and observational bounds. Recent experimental data favour a positive 
cosmological constant, nevertheless since zero cosmological constant is not completely 
ruled out it can be argued that also a tiny negative cosmological constant can be 
acceptable within the experimental uncertainties. The effective cosmological constant 
on the two branes is: 


Am = = -(4.12) 

cosh {kyQ) 

From the above formula we can see that we can make the cosmological constant 
small enough |A 4 d| < if we choose large enough kyo, i.e. kyo > 135. This 

however will make observable deviations from Newtonian gravity at ultra-large scales 
impossible as we will see later. 

To determine the phenomenology of the model we need to know the KK spectrum 
that follows from the dimensional reduction. This is determined by considering the 

linear fluctuations of the metric around the vacuum solution that we found above. We 
can write the metric perturbation in the form; 

y)] dx^dx’' -F dy"^ (4.13) 

where is the four dimensional HdS '4 metric. Here we have ignored the dilaton 
mode that has been calculated in the previous chapter. We expand the held h^y{x,y) 
into graviton and KK plane waves: 


( 4 - 14 ) 

n .=0 


where we demand (V^V^ -|- 2H‘^ — ml) = 0 and additionally = 0 . 
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The function will obey the second order differential equation; 




(4.15) 


After a change of variables and a redehnition of the wavefunction the above equation 
reduces to an ordinary Schrodinger-type equation: 

+ v(z)'^ <i-<">(j) = (4.16) 

where the potential is given by: 






3/c 

T 




cos^ 


tanh{kyQ)5{z) + 


zo) 


sinh(fc(L - yo)) 

cosh(%o) 


5{z - zi) 


(4.17) 


with k dehned as k = —r4—r. The new variables and the redehnition of the wave- 

cosn(fcyo) 

function are related with the old ones by: 


z = sgn(|/)y 
k 


arctan ^tanh( 


k{\y\ - Vol¬ 


ar ctan ^tanh(^^' 


(4.18) 


<i<(">(z) = —i=4.<”l(;,) (4.19) 

Va(2/) 

Thus, in terms of the new coordinates, the branes are at zlq = 0 and zl, 
with the minimum of the potential at zq = | arctan (tanh(^)). Also note that 
with this transformation the point y = oo is mapped to the hnite point Zoo = 
I -I- arctan (tanh(^))]. 

From now on we restrict ourselves to the symmetric conhguration of the two branes 
with respect to the minimum zq (he. the hrst brane at 0 and the second at 2zo ), since 
the important characteristics of the model appear independently of the details of the 
conhguration. Thus, the model has been reduced to a quantum mechanical problem 
with (5-function potentials wells of the same weight and an extra smoothing term in- 
between (due to the AdS geometry). This gives the potential a double “volcano” form 
(see Fig. |4.2|) . An interesting characteristic of this Schrodinger equation is that, as in 
the hat brane case, it has a “supersymmetric” form (see 0) which guarantees that 
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Figure 4.2: The potential V{z) in the ” + +" bigravity model for the symmetric 
conhguration. The (5-function wells correspond to branes. 


it has no tachyonic modes. 

Since the system is compact, i.e. for zl < z^o, it gives rise to a normalizable 
massless zero mode, something that is expected since the volume of the extra dimension 
is hnite. The zero mode wavefunction is given by: 


^(o)(^) = 


A 


lcos(k(zo — 1^1))]^'^^ 


(4.20) 


r-ZL 


where the normalization factor A is determined by the requirement / dz (z) = 1, 

J-zl 

chosen so that we get the standard form of the Fierz-Pauli Lagrangian. 

The form of the zero mode resembles the one of the zero mode of the " H- 

model, i.e. it has a bounce form with the turning at zq (see Fig. [4.3|) . In the case of the 

" H- V" model the cause for this was the presence of the " brane. In the present 

model considering AdS 4 ^ spacetime on the branes we get the same effect. 

In the case that we send the second brane to inhnity {i.e. z —>• Zoo) the zero mode 
fails to be normalizable due to singularity of the wavefunction exactly at that point. 
This can be also seen from (|4.1CI|) which implies that at this limit M* becomes inhnite 
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{i.e. the coupling of the zero mode becomes zero). Thus in this limit the model has no 
zero mode and all gravitational interactions must be produced by the ultralight hrst 
KK mode. The spectrum in this case was discussed by Randall and Karch in . 


Considering the Schrodinger equation for m 7 ^ 0 we can determine the wavefunc- 
tions of the KK tower. It turns out that the differential equation can be brought to a 
hypergeometric form, and hence the general solution is given in terms two hypergeo¬ 
metric functions: 


= cos^/‘^(k(\z\ - zo)) Cl F{an, K, sin^(fc(|z| - Zq))) 

+ C2 I sin(fc(|z| - zo))| F{an + \,bn + |; sin^(fc(|z| - zo))) 


(4.21) 


where 


Cln 



b 


n 



(4.22) 


with the same normalization convention as the zero mode. The boundary conditions 
(he. the jump of the wave function at the points z = 0, zl) result in a 2 x 2 ho¬ 
mogeneous linear system which, in order to have a non-trivial solution, should have 
vanishing determinant. In the symmetric conhguration which we consider, this pro¬ 
cedure can be simplihed by considering even and odd functions with respect to the 
minimum of the potential Zq- 

In more detail, the odd eigenfunctions obeying the boundary condition T*^”)(zo) = 0 
will have Ci = 0 and thus the form: 

= C2CosC‘^{k{\z\ - Zo))|sin(fc(|2;| -Zo))| F{dn + ^,bn + - Zq))) 

(4.23) 

On the other hand, the even eigenfunctions obeying the boundary condition 
'[zq] = 0 will have C 2 = b and thus the form: 

= Cl cos^/^(fc(|z| - Zii))F{an, bn, sm^{k{\z\ - zo))) (4.24) 

The remaining boundary condition is given by: 

3h 

4,(») '(0) + taiih(49„)®<”l(0) = 0 


(4.25) 
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22/0 


o 

o 

o 


m 


Figure 4.3: The graviton (blue line), first (red line) and second (green line) KK states 
wavefunctions in the symmetric " + +" model. All wavefunctions are continuous at 
the position of the bounce. On the right the pattern of the spectrum is sketched. 


and determines the mass spectrum of the KK states. From this quantization condition 

we get that the KK spectrum has a special hrst mode similar to the one of the " H- V" 

bigravity model. For ky^ > 5 the mass of the hrst mode is given by the approximate 
relation: 

mi = 4^3 k (4.26) 

In contrast, the masses of the next levels, if we put together the results for even 
and odd wavefunctions, are given by the formula: 

m„+i = 2^/n{n + 3) k 6“^^° (4.27) 


with n = 1, 2,.... 

We note that the hrst KK state has a diherent scaling law with respect to the 
position of the minimum of the warp factor compared to the rest of the KK tower, 
since it scales as while the rest of the tower scales as Thus, the hrst 

KK state is generally much lighter than the rest of the tower. It is clear that this 

mass spectrum resembles the one of the " H- C' bigravity model. The deeper reason 

for this is again the common form of the warp factor. In both cases the warp factor 
has a minimum due to its bounce form. The graviton wave function follows the form 
of the warp factor, i.e. it is symmetric with respect to zq, while the wavefunction 
of the hrst KK state is antisymmetric in respect to Zq (see Fig.^^BI). The absolute 
values of the two wavefunctions are almost identical in all regions except near zq (the 
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symmetric is nonzero while the antisymmetric is zero at zq). The graviton wavefunction 
is suppressed by the factor cosh“^(%o) at zq which brings its value close to zero 
for reasonable values of kyo. Thus, the mass difference which is determined by the 
wavefunction near zq is expected to be generally very small, a fact which appears as 
the extra suppression factor in the formula of mi in comparison with the rest of 
the KK tower. 

In the case that we consider an asymmetric brane configuration, for example zl > 
2zo, the spectrum is effectively independent of the position of the second brane zl 
(considering kyQ > 5). Thus, even in the case that we place the second brane at z^o, 
i. e. the point which corresponds to infinity in the ^/-coordinates, the spectrum is given 
essentially by the same formulas. This decompactification limit is exactly what has 
been studied in . 


The coupling also of the first KK state to matter on the branes is very different 
from the one of the remaining tower. We can easily calculate it by the same proce¬ 
dure presented in chapter 2 considering the a minimally coupled to gravity matter 
Lagrangian. In the case of scale equivalent positive tension branes, the gravitons will 
couple to matter on them [z = Zbr = 0 or L) with the following couplings: 


— 


jZbr) 

2V2IP 


(4.28) 


In more detail we find that the zero mode and the first ultralight KK state couple 
with the same strength on the positive tension branes: 


Gq — Gi — 


2V2M, 


(4.29) 


while the coupling of the remaining tower is exponentially suppressed (it scales as 
e~^y°). The reason for the strong coupling of the first KK state in comparison with 
the remaining tower is readily understood because, being dominantly a bound state 
of the volcano potential on the positive tension brane, it is largely localized on it. 

Exploiting the different mass scaling of the first KK relative to the rest we can 

ask whether it is possible to realize a bigravity scenario similar to that in ” H- \-" 

bigravity model. In that model by appropriately choosing the position of the negative 
tension brane, it was possible to make the first KK state have mass such that the 
corresponding wavelength is of the order of the cosmological scale that gravity has 
been tested and at the same time have the rest of the KK tower wavelengths below 
0.2mm. In this scenario the gravitational interactions are due to the net effect of the 
massless graviton and the first ultralight KK state. A similar situation can be realized 
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in the " + +" model as well. From ( |4.29|) it can be understood that in the symmetric 
" + +" conhguration the massless graviton and the special KK state contribute by the 
same amount to the gravitational interactions. In other words: 


1 _ 1 1 


Mpi = 




(4.30) 


However, if our aim is to see modihcations of gravity at ultralarge distances, this 
is impossible because the Compton wavelength of our ultralight graviton will be 
times bigger than the horizon H~^ of the AdS 4 ^ space on our brane due to (|4.12|) , 
(|4.26|) . The “Hubble” parameter follows m 2 rather than mi. Moreover, in the present 
model, the fact that the effective four dimensional cosmological constant should be set 
very close to zero, requires that the warp factor is constrained by kyo >135 and thus, 
in this case, the spectrum of the remaining KK states will be very dense (tending to 
continuum) bringing more states close to zero mass. The KK states that have masses 
which correspond to wavelengths larger than 1mm have sufficiently small coupling so 
that there is no conflict with phenomenology (the situation is exactly similar to the 
RS2 case where the coupling of the KK states is proportional to their mass and thus it 
is decreasing for lighter KK states). The fact that the spectrum tends to a continuum 
shadows the special role of the hrst KK state. 

It is interesting to note that at the limit where the minimum of the warp factor is 
sent to inhnity (zoc) the special behaviour of the hrst KK persists and does not catch 

the other levels (by changing its scaling law) as was the case in " H- \-" model of 

chapter 2. This means that the limit z ^ z^o will indeed be identical to two RS2, 
but on the other hand we should emphasize that what we call graviton in the RS2 
limit is actually the combination of a massless graviton and the “massless” limit of the 
special hrst KK state. This “massless” limit exists as we will see in the next section 
and ensures that locality is respected by the model, since physics on the brane does 
not get ahected from the boundary condition at inhnity. 

Finally, let us note that one can build a crystalline model as in the hat brane 
case but without the need to introduce negative tension branes. In this case one 
expects again to get a band structure for the KK continuum with the width of the 
hrst band (starting from zero) much smaller than the width of the hrst forbidden 
zone. However, as in the bigravity case the width of the hrst band which will mediate 
normal four dimensional gravity at intermediate distances will be linked to the remnant 
cosmological constant on the branes and quite probably one will get modihcations of 
gravity for distances much bigger than the AdS horizon. The analogue, however, of 
the quasi-localized model cannot be realized in the AdS brane case as it can be seen 











2 Graviton propagator in dS 4 and AdS^ space 


67 


from (|4.4|), ([4.51) since an asymptotically constant warp factor is forbidden. 


2 Graviton propagator in dSi and AdSi space 


In this section we will demonstrate how the van Dam-Veltman-Zakharov no-go theorem 
for the non-decoupling of the extra polarization states of the massive graviton can be 
evaded in dS '4 or AdS 4 space. For the case of dS '4 there has been a proof in but 

we will present an alternative one in this section together with the case of AdS 4 for 
the sake of completeness. In more detail, as we have discussed at the end of chapter 
2 , in flat spacetime the m —»■ 0 limit of the massive graviton propagator does not give 
the massless one due to the non-decoupling of the additional longitudinal components. 
This generates the well known discontinuity between massive and massless states. 
Considering the massive graviton propagator in (iS '4 or AdS 4 spacetime we can show 
that this result persists if m/H —> 00 where H is the “Hubble” parameter, i.e. the 
discontinuity is still present in the m —> 0 limit if it happens that m/H 00 . However, 
in the case that m/H —> 0, we will explicitly show that the m —»■ 0 limit is smooth. 
This is an important result since it gives us the possibility to circumvent the van 
Dam - Veltman - Zakharov no go theorem about the non-decoupling of the extra 
graviton polarizations. Thus, in the limit that m/H —> 0 all the extra polarizations of 
the graviton decouple, giving an effective theory of a massless graviton with just two 
polarization states. 

Our aim is to hnd the form of the massless and massive graviton propagators in the 
case of (iS '4 and AdS 4 spacetime with arbitrary “Hubble” parameter H and graviton 
mass m. We wish to examine the behaviour of these propagators in the limit where 
H > 0 and the limit where the mass of the graviton tends to zero. For simplicity we 
will do our calculations in Euclidean (iS '4 or AdS 4 space. We can use for metric the 
one of the stereographic projection of the sphere or the hyperboloid[J 


ds^ = 


5 


/il/ 




dx^dx'' = g^^j^dx^dx’' 


(4.31) 


where x"^ = d^j^^x^x'^ and the scalar curvature is i? = 4=12 H^. From now on, the 
upper sign corresponds to AdS 4 space while the lower to (iS '4 space. The fundamental 
invariant in these spaces is the geodesic distance p^{x,y) between two points x and 
y. For convenience, we will introduce another invariant u which is related with the 
geodesic distance by the relation u = cosh(i7/i) — 1 for AdS 4 {u G [ 0 , cx))) and the 

^Note that and whose results we use in the following have a different metric convention, 
but this makes no difference for our calculations. 
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relation u = cos{Hfi) — 1 for dS 4 {u G [—2, 0]). In the small distance limit u 
This background metric is taken by the the Einstein-Hilbert action; 




Sg= / d^x^{2MlR-K) 


(4.32) 


where the cosmological constant is A = =pl2 and M 4 the four dimensional 

fundamental scale. The spin-2 massless graviton held can be obtained by the linear 
metric huctuations ds"^ = dx^dx’^. This procedure gives us the analogue of 

the Pauli-Fierz graviton action in curved space: 


^0 

2M| 





K 


± -H^ 
2 


:h^ + Kuh>^ 


(4.33) 


The above action is invariant under the gauge transformation 

which guarantees that the graviton has only two physical degrees of freedom. This is 
precisely the dehnition of masslessness in dS '4 or AdS^ space (for example see [|^ and 
references therein). 

The propagator of the above spin-2 massless held can be written in the form: 




+ d^dy>udydf,>u)G'^{u) + gf,ugyv>E^{u) + D[-■ ■] 


(4.34) 


where dy = The last term, denoted £)[■■■], is a total derivative 

and drops out of the calculation when integrated with a conserved energy momentum 
tensor. Thus, all physical information is encoded in the hrst two terms. 

The process of hnding the functions and is quite complicated and is the result 
of solving a system of six coupled diherential equations . We will present here only 
the diherential equation that G^ satishes to show the diherence between AdS 4 ^ and (iS '4 
space. This equation results from various integrations and has the general form: 


u{u + 2)G'°(m)" + 4(m + 1)G'°(m)' = Cl + G 2 U (4.35) 


where the constants Ci and C 2 are to be hxed by the boundary conditions. For the 
case of the AdS 4 ^ space [^], these constants were set to zero so that the function 
vanishes at the boundary at inhnity {u —*■ cxd). Using the same condition also for the 






2 Graviton propagator in dS 4 and AdS^ space 


69 


function, the exact form of them was found to be: 


E^(u) 


2(u + 1) 

[u{u + 2) 

r2(u + l) 


log 


M + 2 


u 


Svr^ [u{u + 2) 


+ 4(m + 1) - 2(m + 1) log 


u + 2 


u 


(4.36) 


For the case of the dS^ space we iterated the procedure of imposing the condi¬ 
tion 1^ that the and E^ functions should be non-singular at the antipodal point 
(m = —2). The constants Ci and C 2 were kept non-zero and played a crucial role in 
hnding a consistent solution. It is straightforward to hnd the full expression of these 
functions, but we only need to know their short distance behaviour. Then with this 
accuracy the answer is: 


G^u) 

E\u) 


1 




u 


+ log(-M) 


+ ■ ■ ■ 




- + 2{u + 1)^ log(-M) 
u 


+ 


(4.37) 


If we dehne 11° (u) = then for short distances -C 1) where 

we get: 


u 




d^d^iu 

TH^d^y, 

G\u) 

1 


n°(u) 

-1 


(4.38) 


and so we recover the short distance limit of the massless flat Euclidean space propa¬ 
gator: 

^~ (4.39) 

Of course this is just as expected. 

In order to describe a spin-2 massive held it is necessary to add to the above action 
a Pauli-Fierz mass term: 



2 M| 


So m? 

2Ml ~ ~T 




(4.40) 


By adding this term we immediately lose the gauge invariance associated with the 
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dS 4 or AdS 4 symmetry group and the graviton acquires five degrees of freedom and 
becomes massive. 

The propagator of this massive spin-2 field can again be written in the form: 


y) = {df,dpud^d^>u -h d^,d^mdydpu)G'^{u) + g^u9yu'E'^{u) + £>[•■■] (4.41) 


The last term of the propagator in (|4.41|) , denoted /![■■■], is again a total derivative 
and thus drops out of the calculation when integrated with a conserved T^y. 

At his point we should emphasize that in case of an arbitrary massive spin -2 field, 
the absence of gauge invariance means that there is no guarantee that the field will 
couple to a conserved current. However, in the context of a higher dimensional theory 
whose symmetry group is spontaneously broken by some choice of vacuum metric, the 
massive spin-2 graviton KK states couple to a conserved T^y. One can understand 
this by the following example. Consider the case of the most simple KK theory, the 
one with one compact extra dimension. By the time we choose a vacuum metric e.g. 
9mn — diag 1), the higher dimensional symmetry is broken. If we denote the 
graviton fluctuations around the background metric by h^y, and ^ 55 , there is still 
the gauge freedom: 


5h^y = d^iy -F dyi^ 

= dfj,^5 + d^^n 

= 2d^^^ 


(4.42) 


If we Fourier decompose these fields, their n-th Fourier mode acquires a mass oc 
n with n = 0,1, 2,..., but there is mixing between them. This means for example 

that is not a massive spin-2 eigenstate etc.. However, we can exploit the gauge 
transformations (|4.42|) to gauge away the massive and and construct a pure 
spin-2 field (see for example |^] and references therein). For a comprehensive account 

of KK theories see . The new massive spin-2 field p^^y is invariant under ([4.42|) and 
so its Lagrangian does not exhibit a gauge invariance of the form 5p^j_y = d^j_Xv + dyXii- 
However, since is originates from a Lagrangian that has the gauge invariance ([4.42 ), 
it is bound to couple to a conserved T^y. The argument goes on for more complicated 
choices of vacuum metric as for example warped metrics. 

Again the functions G™ and result from a complicated system of differential 


equations |69]. In that case, the differential equation that G™' satisfies is: 


u{u + 2)G^{u)'' + 4{u + l)G”^(n)' ^ G™(n) = Gi + Gsw 


(4.43) 
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where the constants Ci and C 2 are to be hxed by the boundary conditions. For the 
case of the AdS 4 space [^, these constants were set to zero so that the function 
vanishes at the boundary at inhnity. Imposing additionally the condition of fastest 
falloff at inhnity (m —>■ 00 ) Q. the exact form of the G™' and function was found 
to be: 


G^(u) = 
E’^(u) = 


X 


r(A)r(A - 1 ) 
167r2r(2A - 2)Fr2 



A 

F(A,A 


l,2A-2,--) 

u 


2 V{A-1)H^ 

3 167r2r(2A-2)[2 +(m/iJ)2] ^ 


3[2 + {m/Hf]r{A - 2)u‘^F{A - 1, A - 2,2A 
< -3(M + l)MF(A-l,A-l,2A-2,-f) 

+ [3 + {m/Hf]T{A)F{A, A - 1, 2A - 2, -^) 



(4.44) 


where A = | + ^ . 

For the case of the dS 4 space we iterated the procedure of imposing the condi¬ 
tion that the and functions should be non-singular at the antipodal point 
{u = —2) and also hnite as m —> 0. Again we kept the constants C 4 and C 2 non-zero 
to obtain a consistent solution. It is straightforward to hnd the full expression of these 
functions, but we only need to know their short distance behaviour. Then with this 
accuracy the answer is: 


G^{u) 

E^{u) 


X 


F(A)F(3-A) 


167r2iy2 

2 F(A)F(3 - A)Ff2 

'3 47r2[2 - [m/Hf] 


7/ j- 2 

F(A,3-A,2,^^)-l 


+ ■ ■■ 


-3[2 - (m/i7)2] 
-3 (m + 1) 


X 


2(u+2) 


(A-i)(a^^(^-1>2-A,2,h^)+ , 

2 -^F(A-1,2-A,1,^) + (m + 1)' 


_(A-l)(A-2) 

+ [3 - [m/Hf] [F(A, 3 - A, 2, ^) - l] 


+ 


where A = | -|- ^ ^y9 — 4:{m/H)^. 

If we dehne 11"* (m) = then for short distances <C 1) where m —> 0 

we get: 




1 


(4.45) 
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n" 


u 


2 

3 


3±(A 


It is interesting to consider two massless flat limits. In the first one m 
77 —> 0 while m/H 


(4.46) 


0 and 


cxo. In this case, from ([4.46|) we see that we recover the 
Enclidean propagator for a massive graviton in flat space: 


1 2 


(4.47) 


This is in agreement with the van Dam - Veltman - Zakharov theorem. The second 
limit has m ^ 0 and 77 —> 0 bnt m/H ^ 0. In this case the propagator passes 


smoothly to the one of the flat massless case ( [4.39| ): 




' 4 ^ 2^2 




(4.48) 


This is in contrary to the van Dam - Veltman - Zakharov discontinnity in flat space. 

In general, we may consider the limit with m/H finite. Then for small m/H the 
contribntion to the strncture is —1 ± m?/QH'^. Since observations agree to 

0.02% accuracy with the prediction of Einstein gravitational theory for the bending 
of light by the sun, we obtain the limit ^ < 0.04. However, we should point out that 
the above results are meaningful only for the AdS case because in the dS there are no 


unitary spin-2 representations in the region 0<m/77<2asit was shown in . 


3 Discussion 

Let us now go back to the ” + +" bigravity model and see how the extra polarization 
states problem is evaded in the absence of a ghost radion mode. In this model we have 
a clear realization of the case where m/H —> 0. Indeed the massive graviton, which 
in the symmetric configuration is responsible for half of the gravitational interactions 
that we feel in intermediate distances, has this characteristic ratio to be mi/H ~ 

This easily accommodates the bound that we presented earlier. Then, the Euclidean 
propagator (in configuration space) of the massive KK states for relatively large yn 
will be given by: 




dvr^/i- 


:(V<^-' + W - (l - (4.49) 
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where ju is the geodesic distance between two points. In the above, we have omitted 
terms that do not contribute when integrated with a conserved For ky^ > 2.3 
there is no problem with the bending of light since the extra polarization states are 
practically decoupled. Additionally, as we discussed earlier the large scale modihca- 
tions of gravity are impossible to be observed because the Compton wavelength of our 
ultralight graviton is exponentially larger than the horizon H~^ of the AdS^ space on 

our brane. The facts that there is no observable deviation from normal four dimen¬ 
sional Einsteinian gravity in this model and that there is a cosmological constant of 

negative sign present, make this model phenomenologically disfavoured. 

Let us now make some further comments regarding the (dis)continuity of the gravi¬ 
ton propagator. It was shown in |]71| , |75| that the extra polarizations are resurrected in 
quantum loops and thus the discontinuity reappears in the quantum level. However, 
purely four dimensional theory with massive gravitons where the mass of the graviton 
is added by hand is not well-dehned as it was shown in . If instead one starts from a 
higher dimensional theory in which the mass of the gravitons is generated dynamically, 
as we are doing in the models we have presented so far, and if one takes into account 
the whole KK graviton tower in the quantum loops, this quantum discontinuity is 
expected to be absent. It is worth noting that in [j^ it was suggested that the extra 
polarization states of the massive graviton in the one AdS^^ brane model [J3| can be 
gauged away and thus the discontinuity is absent. 

Moreover, one expects that the smoothness of the limit m —*■ 0 is not only a prop¬ 
erty of the AdS '4 space but holds for any background where the characteristic curvature 
invariants are non-zero. For physical processes taking place in some region of a curved 
space with some non-zero characteristic average curvature, the effect of graviton mass 
is controlled by positive powers of the ratio m?/B? where is a characteristic curva¬ 
ture invariant (made from Riemann and Ricci tensors or scalar curvature). 

Finally, the conjecture that there is a smooth limit for phenomenologically observ¬ 
able amplitudes in brane gravity with ultralight gravitons is also supported by results 


based on a very interesting paper [^. In that paper it was shown that the argument 
of van Dam, Veltman and Zakharov that no mass for the graviton is allowed in flat 
background is not entirely robust. In more detail, is was shown that there is a smooth 
limit for a metric around a spherically symmetric source with a mass M in a theory 
with massive graviton with mass m for distances smaller than a critical one which 
is m~^{mM/. This critical distance for a graviton of wavelength of the size 
of the observable universe and for a massive body of the mass of the Sun, is of the 
order of the solar system size. In order to see this smooth limit one should perform 
a non-perturbative analysis since in this regime the lowest order approximation 
breaks down. The discontinuity reappears for distances larger than this critical dis- 
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tance where the linear approximation is valid. This breakdown of perturbativity, is 
due to the fact that the lowest order approximation calculation does not take into 
account the characteristic mass scale of the problem which in this case is the mass of 
the Sun. Since the bending of the light is observed in the hrst non-perturbative region, 
according to []^ it is impossible to rule out the massive graviton proposal. 

flowever, the non-perturbative solution discussed in p0| , [75[] was found in a limited 
range of distance from the center and it is still unclear if it can be smoothly continued 
to spatial inhnity (this problem was stressed in [^). Existence of this smooth contin¬ 
uation depends on the full nonlinear structure of the theory. If one adds a mass term 
by hand the smooth asymptotic at inhnity may not exit. However, it seems plausible 
that in all cases when modihcation of gravity at large distances comes from consistent 
higher-dimensional models, the global smooth solution can exist because in this case 
there is a unique non-linear structure related to the mass term which is dictated by 
the underlying higher-dimensional theory. 







Chapter 5 

Multigravity in six dimensions 


As we have discussed, in five dimensions it is impossible to have fiat brane multigravity 
models without negative tension branes. For this reason, we will consider models in six 
dimensions. The literature on six dimensional constructions is already quite rich [|7^|80| , 
^,|8^ , ^,|8^ , ^,|86| ,^ , |88|,^ , |^,^ , P^,[^ , |9l|,[^ , |^ [^. In this chapter we explicitly show 

that the difficulties of the five dimensional models can be evaded in six dimensions. 
It is possible to construct multigravity models with only fiat positive tension branes. 

This is due to the fact that in six dimensions the internal space is not trivial and 
its curvature allows bounce configurations with the above feature. The branes which 
localize gravity in this setup are four-branes, but one of their dimensions is compact, 
unwarped and of Planck length. Thus, in the low energy limit the spacetime on the 
brane appears three dimensional. In order that these constructions are realized it is 
crucial that the tensions and/or the bulk cosmological constant are anisotropic. The 
five dimensional constructions in this setup come into two types. One of them involves 
conical singularities at finite distance from the four branes. These conical singularities 
support three-branes which can be of positive tension if one has an angle deficit, of zero 
tension if one has no angle deficit and of negative tension if one has an angle excess. The 
other type of the multigravity models has no conical singularity at all. It is crucial to 
note that there is no five dimensional effective theory for these constructions, otherwise 
one would get all the problems faced in the five dimensional constructions. The low 
energy effective theory is directly four dimensional. These constructions give for the 
first time a theoretically and phenomenologically viable realization of multigravity. 

In the following section we present the cylindrically symmetric single four-brane 
warped models and show the cases where gravity is localized on the four-branes. In 
section two we paste two of the single four-brane solutions and obtain the double 
four-brane bigravity models. In section three we generalize the quasi-localized and 
crystalline constructions in six dimensions. Finally, we discuss the weaker energy 
condition in these six dimensional configurations. 
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a, 7 



*■--4 



P 


Figure 5.1: The minimal single four-brane model warp factors cr(p) and y{p)- 


1 Single brane models in six dimensions 

At first we will discuss the single brane solutions in six dimensions to get an insight 
on the more complicated multigravity conhgurations. In the following subsection we 
will review the minimal model where the bulk cosmological constant is isotropic and 
then a more generalized model with arbitrary cosmological constant components along 
the four dimensional and the extra dimensional directions. In the following we will 
implicitly assume orbifolding around the four-brane. 

1.1 The minimal single brane model 

The simplest single brane model consists of a four-brane embedded in six dimensional 
AdS space [^. One of the longitudinal dimensions of the four brane is compactihed 
to a Planck length radius R while the dimensions transverse to the four-brane com¬ 
pact dimension are inhnite. We consider the following ansatz where four dimensional 
Poincare invariance is respected^: 


ds"^ = (j{p)rj^j_ydx^dx^ + dp^ + '-y{p)d9‘^ 


(5.1) 


where 6 is the compactihed dimension with range [0, 27r] and p is the inhnite radial 
dimension. 

It is straightforward to solve the Einstein equations for the bulk energy momentum 


tensor = — with the brane eontribution: 



(5.2) 


^Note that in this thesis we use different metric signature and different definition of the funda¬ 
mental scale from [^. 
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From the form of the Einstein eqnations given in Appendix C we see that the 
solution for the two warp factors is: 


a(p) = e 


— p-kp 


7(p) = 


(5.3) 


with /c^ = — yq^, where M is the six dimensional fundamental scale and the arbitrary 
integration constant R is just the radius of the four-brane (see Fig. ^.ll) . The Einstein 
equations require the usual hne tuning between the bulk cosmological constant and 
the tension of the four brane: 

8 A 


V = - 


5k 


(5.4) 


Let us note at this point that in this hne tuning was absent because a smooth local 
defect was considered instead of a four-brane. In this case, the hne tuning emerges 
between the diherent components of the defect energy momentum tensor. The physics 

of the four-brane idealization and the one of the defect model is the same. 

The four dimensional KK decomposition can be carried out as usual by considering 
the following graviton perturbations: 


= a{p) + h^u{p, 6, x)] dxPdx'" + dp^ + p{p)d6" 


(5.5) 


Here, as well as throughout this chapter, we have ignored the modulus associated 

to the radius R of the four-brane. This will be a massless scalar and in order not 
to give rise to an additional long range force it should be given mass through some 
stabilization mechanism. 

We expand the graviton perturbations in a complete set of radial eigenfunctions 
and Fourier angular modes: 


h^^{p,e,x) = ^(j)(n,i){p)e'^^h\llfHx} 




(5.6) 


n,l 


The diherential equation for the radial wavefunctions (j) is: 


+ (r “ ^ 1 ^ 0 


(5.7) 


with normalization dpa^cjy^cjyn = 5mn- We can convert this equation to a two 
dimensional Schrodinger-like equation by the following redehnitions: 




(5.8) 
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Figure 5.2: The form of the effective potential 14// for different angular quantum 
numbers I for the minimal single four-brane model. The h-function in the origin is 
omitted. 


so that 


1 / ■^ \ 777 ^ '' 21 


3k 


2i?2 


S(z) (5.9) 


From the form of the potential (see Fig. |5.2|) we can easily deduce that the angular 
excitations spectrum will consist of continunm starting from a gap of the order and 
thus can be safely ignored. For the s-wave (/ = 0) there is a normalizable zero mode 
which is a constant in the p coordinate, i.e. 0o = const. The KK tower for the s-waves 
will again form a continunm but this time gapless with wavefunctions given by 




f2m\ ( 2m k 

(x) (x" 


2m k 


(2m\ k 


(5.10) 


The correction to Newton’s law due to the s-modes can be easily calculated and 
one hnds that it is more suppressed than the RS2 model In six dimensions 

the correction reads: 


Al/ = 



(5.11) 


1.2 The generalized single brane model 

We can now try to hnd more general single four-brane solutions by relaxing the re¬ 
quirement of isotropy of the bnlk and brane energy momentum tensors. In contrast 
to the hve dimensional case, this is possible in six dimensions because there are two 
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independent functions in the metric (see Appendix C). We will consider the following 
anisotropic bulk energy-momentum tensor: 


^(B) N 
M 


An 5 


U 



and allow for an anisotropic brane tension of the form: 


rp{br) N 

^ M 


-i(p) 



(5.12) 


(5.13) 


This anisotropy can be due to different contributions to the Casimir energy in the 

(see also |^) or due to a background three-form gauge 


different directions 


held with non-zero held strength 


If we dehne the parameter a = then the one four-brane solutions for the warp 
factors are: 

a(p) = e-'^P , 7(p) = /?2g-|(5a-6)p 



Figure 5.3: The generalized single four-brane warp factors cr(p) and 7 (p) for diherent 
values of the parameter a. 
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Figure 5.4: The form of the effective potential Vg/f for a > 2 for different angular 
quantum numbers 1. The ^-function in the origin is omitted. 


with = — ^^ 4 , where again R is the radius of the four-brane (see Fig. b.3|) . The 
flatness of the brane is achieved by the following hne tunings of the the parameters of 
the bulk and brane energy momentum tensors: 


5k ’ 16 


A + 12 ^ 


(5.15) 


Let us note that the two components of the brane tension have the same sign as 
long as a > — |. For « > | the internal space is shrinking as in |^, for « < | the 


internal space is growing and for a = | the internal space is a cylinder (see Fig.^ 
Proceeding with the usual KK decomposition we hnd the following equation for 



Figure 5.5: The form of the effective potential 14// for a < 2 for different angular 
quantum numbers /. The (5-function in the origin is omitted. 
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the radial wavefunction 0: 


5a+ 10 
8 


k<^' + 



) 


e’^Pcj) = 0 


(5.16) 


with normalization dpay/y(l)m4>n = dmn- With the wavefunction redefinition and 
the coordinate change discussed in the previous subsection, we obtain the following 
effective potential for the two dimensional Schrodinger-like equation: 


V.„(z) 


3(5a + 4)/c^ 

64(t + 1)' 


P Pkz 

2 ^ W 




(5.17) 


For the s-wave there is a normalizable zero mode when a > — | which is a constant 

5 

in the p coordinate, i.e. 0o = const. The KK tower for the s-waves will again form a 
continuum and their wavefunctions are given by: 


0m 







2m 

2 ^ 


(5.18) 


where u = |(a + 2). The correction to the Newton’s law due to the s-modes can be 
easily calculated and one finds that it is negligible as long as a ^ 


Al/ 



(5.19) 


As far as the angular excitations I ^ 0 are concerned, in the case where a > 2 the 
potential diverges at inhnity (see Fig. |^) , so the spectrum will be discrete starting 

from a scale bigger than ^ and thus can be safely ignored. On the other hand, 
the potential for the case a < 2 asymptotically vanishes (see Fig. ED and so the 
spectrum will be gapless continuum. One should examine the contribution of the 
angular excitations to the Newtonian potential in this case. 


2 Bigravity in six dimensions 


The flat single four-brane models considered in the previous sections have the charac¬ 
teristic that gravity is localized on the brane in the same way as in the five dimensional 
analogue [1^. In this section we will show how we can construct realistic multi¬ 
localization scenarios for gravity by sewing two single four-brane solutions together. 
We will present two explicit examples of compact double four-brane models where bi- 
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gravity is realized. Firstly, we will consider a double four-brane bigravity model which 
in general contains an additional three-brane associated with the existence of a conical 
singularity. This can be done by allowing for an anisotropic four-brane tension. In the 
following subsection we demonstrate how one can avoid the conical singularity by also 
allowing for anisotropic bulk cosmological constant. In the following we will consider 
the orbifold geometry jZ^ in which the two four-branes lie on the hxed points. 


2.1 The conifold model 

We are interested in a double four-brane generalization of the minimal single four-brane 
model. In order to achieve this, we still assume an isotropic bulk energy momentum 

tensor but we do not impose any constraints on the four-brane tension: 


rpibr) N 

^ M 


-6{p) 



(5.20) 


The Einstein equations for u{p) and y{p) in this case give the following solutions for 
the warp factors: 


_ cosh^/^ [\k{p - po)] _ [ffcpo] sinh^ [\k{p - Po)] g.x 

cosh^/® [|/cpo] ’ ^ sinh^ [|/cpo] cosh®/*^ [|/c(p - po)] 


with where we have normalized (t( 0) = 1 and 7 ( 0 ) = (see Fig. [^) . 

From the above relations it is obvious that both u{p) and y{p) have a bounce form. 
However, we note that 7 (po) = 0, that is, y{p) is vanishing at the minimum of the 
warp factor ct(p). This is a general characteristic of the solutions even when the branes 


are non-flat. From ( |C.3|) and (|(J.4|) we can easily show that y{p) = (where C 

is an integration constant) which implies that whenever we have a bounce in the warp 
factor the function y{p) will develop a zero. 

In order to examine the nature of this singularity we study the form of the metric 
at the vicinity of the point p = po along the lines of ||100| , p.01|] . Taking in account that 
in this limit we have a{p) 6 ^ and 7 (p) —*• (3‘^{p — po)^ the metric becomes: 


ds"^ = 


p^udx^dx’' + dp^ + j3‘^{p — po)^d6‘^ 


where 6 ^ 


cosh "‘Z® [f fcpo] 


(5.22) 
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Figure 5.6: The conifold bigravity model warp factors a(p) and j(p) for the symmetric 
conhguration where L = 2po- Both the a(p) and 7 (p) have a bounce at po- However, 
jIp) vanishes at this point. This point in general corresponds to a conical singularity. 

From the form of the metric, it is clear that for general values of the jS parameter 
there will be a conical singularity with a corresponding dehcit angle 5 = 27r(l — /?). 
The existence of this conifold singularity is connected to the presence of a 3-brane at 
p = Pq. In order to hnd its tension one has to carefully examine the Einstein tensor 
at the vicinity of the singularity. For this reason we write the metric for the internal 
manifold in the conformally flat form: 

= b^rj^ydx^dx’^ + f{r){dr'^ + r‘^d6‘^) (5.23) 

with /(r) = and p — po = . In these coordinates it is easy to see how the 

three brane appears at the conifold point. The Einstein tensor can be calculated for 
P Po- 

VHog(/(r)) jO \ 

2/(r) ^ 'Irv \ 

0 / 

where is the flat two dimensional Laplacian. Now, given that log(r) = 


(5.24) 


Rmn — -GmnR — 
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and by comparing with: 


R 


MN 


^CmnR ^^4 




(5.25) 


where ^ (see Appendix of [0), we hnd that the tension of the 3-brane is: 


i/3 = 4(1 - /5)M^ = -—6 

71 


(5.26) 


Thus, if there is angle dehcit <5 > 0 (/? < 1) the tension of the brane is positive, 
whereas if there is angle excess 5 < 0 (/3 > 1) the tension of the brane is negative. At 
the critical value /3 = 0 there is no conical singularity at all and we have a situation 
where two locally flat spaces touch each other at one point. 

In the following we will concentrate on the symmetric case, that is, the four-branes 
will be considered placed at symmetric points with respect to the minimum of the 
warp factor po, i.e. L = 2po- For the previous solution to be consistent the tensions 
of the four-branes for the symmetric conhguration must satisfy: 




8 A 

- tanh 
5k 




8A2 1 
5^Ve 


(5.27) 


The brane tensions for both branes are identical. Thus, the above construction 
consists of two positive tension four-branes placed at the end of the compact space 
and an intermediate three-brane due to the conifold singularity with tension depending 
on the parameters of the model. In the limit po 00 we correctly obtain two identical 
minimal single four-brane models for the case where 6 —>■ 277 (/3 —> 0). 

The differential equation for the radial wavefunction 0 of the graviton excitations 


is: 


+ 2 


cr 7 

-1" '7~ 

a 47 






(5.28) 


with normalization dpa^/y(()rn(t>n = dmn- 

There is an obvious normalizable zero mode with 00 = const, and a tower of 
discrete KK states. It is easy to demonstrate that for 5 > 0 (0 < 1) this model has 
an ultralight state, leading to a theory of bigravity without negative tension branes. 
This follows considering the limit po ^ 00 when the theory describes two identical 
single four-brane models each of which has a massless graviton. For hnite po only 
their symmetric combination remains massless but their antisymmetric combination 
acquires an anomalously small mass, because the difference in the modulus of the 
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Figure 5.7: The form of the effective potential 14// for different angular quantum 
numbers I for the conifold bigravity model in the symmetric conhguration. The 6- 
functions in the positions of the four-branes are omitted. 


wavefunctions between the massive and the massless states is signihcant only near 
Po where the wavefunctions are exponentially smal|^ (see ||102|| for discussion of this 
phenomenon for helds of all spin). 

We can transform this differential equation to a two dimensional Schrodinger-like 
one by the coordinate change ^ = g{z) and the usual wavefunction redehni- 

tion. Then the effective potential reads: 


yeff = 


15 


dzQ 

9 


^ dlg _ 3 d^gd^y f 

A g 8 gy 2yg^ 


(5.29) 


We cannot write an explicit analytic formula for of the above potential in the z- 
coordinates because the coordinate transformation is not invertible analytically. How¬ 
ever, since the transformation is monotonic, we can easily sketch the form of the 
potential (see Fig. |5.7|) by calculating it in the p coordinates. From this procedure we 
see that the potential for the s-wave is hnite at the conical singularity po but has a 
divergence for all the angular excitations with I ^ 0. This means that only the s-wave 
excitations with communicate the two parts of the conifold and the other excitations 
will be conhned in the two semicones. The spectrum of KK states in this case will be 

^ There is, however, a mismatch of the degrees of freedom of the massive and massless graviton, 
which may not be phenomenologically dangerous, but still raises the issue of locality in the infinite 
four-brane separation limit. 













Chapter 5: Multigravity in six dimensions 


discrete. 

The construction may be readily modified to obtain a six-dimensional analogue of 


the locally localized model This is done by considering the asymmetric situation 
in which one of the fixed point four-branes is moved to infinity (in the asymmetric 
case the tensions of the two four-branes are different). Then gravity on the four-brane 
at p = 0 will be mediated by only the ultralight state since the graviton zero mode 
will not be normalizable. 


2.2 The non-singular model 


The appearance of a conifold singularity can be avoided if instead of using the minimal 
single four-brane solution to build the double four-brane solution, one uses the gener¬ 
alized single four-brane model. By allowing for an isotropic bulk cosmological constant 
one can arrange that the function 7 (p) does not develop a zero. In this case the (0, 9) 
and (p, p) part of Einstein equations in the bulk, (|U.3| ) and (|C.4| ), have solutions: 


ct(p) = 


cosh^/^ [iHp- po)] 

cosh^/® [|/cpo] 


7(p) = A 


sinh° [\k{p - po)] 
cosh®/'^ [\k{p - Po)] 


(5.30) 


where A is a constant and a = The (p, z/) component of the Einstein equations 




(|C.5|) , however, restricts the possible values of a to a = 2 or cc = 0. The case of a = 2 
corresponds to the previous conifold model. For the case a = 0, which corresponds to 
the choice Ap = 0, the solution is: 




cosh°''^ [^kpo] 
cosh»/^ Yik(p - p„)] 


(5.31) 


with V — — , where we have normalized c 7 ( 0 ) — 1 and 7 ( 0 ) — (see Fig. |^). 

From the above the absence of any singularity is obvious since y{p) does not vanish 
at any finite value of p. 

The four branes will appear as ^-function singularities at the points where we will 
cut the space in the p direction. In the following we will concentrate on the symmetric 
case, that is, the branes will be considered placed in symmetric points in respect to 
the extremum of the warp factor po, i.e. L = 2po. In order for the above configuration 
to be realized, the components of the four-brane tensions and the bulk cosmological 
constant must be tuned to give: 




-^tanh 

5k 



Vo = 



Ao — 



(5.32) 
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Figure 5.8: The non-singular bigravity model warp factors a(p) and 7 (p) for the 
symmetric conhguration where L = 2po- The a(p) warp factor has exactly the same 
form as in the conifold model. The 7 (p) warp factor however now has an inverse 
bounce form. In this case j(p) does not vanish anywhere and thus the model is free 
of singularities. 


The brane tensions for both branes are identical. Thus, we have constructed a 
compact model with two positive tension flat branes with a bounce in the warp factors. 
Note the extra condition relating the different components of the bulk cosmological 
constant; this anisotropy in the bulk cosmological constant is what allows us to have a 
non-singular solution. In the limit po ^ oc we obtain two identical generalized single 
four-brane a = 0 models. 

The differential equation for the radial wavefunction 0 of the graviton excitations is 
given by (|5.28|) for the relevant a(p) and 7 (p) functions. Although it is difficult to hnd 
an analytical solution of the corresponding differential equation for the KK states, the 
existence of a light special KK is assured from the locality arguments discussed above. 
Since at the inhnite brane separation limit we recover two identical generalized single 
four-brane models (for a = 0) where each of them supports a massless zero mode, at a 
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Figure 5.9: The form of the effective potential 14// for different angular quantum 
numbers I for the non-singular bigravity model in the symmetric conhguration. The 
h-functions in the positions of the four-branes are omitted. 

hnite separation conhguration the one of the zero modes will become the special light 
KK state. 

By making an appropriate change of coordinates and a redehnition of the wavefunc- 
tion we can bring the previous differential equation in a two dimensional Schrodinger 
form with the effective potential (|5.29|) that is plotted in Fig. ^.9| . Note that in this 
case the angnlar excitations feel the whole region between the two positive branes and 
thus communicate physics between them. The spectrum of KK states in this case will 
again be discrete. 

3 Multigravity in six dimensions 

In this section we will show how we can obtain the analogues of the GRS model of 
quasi-localized gravity and of the crystal universe in six dimensions, in both cases 
without the need of introducing moving negative tension branes. 

3.1 Quasi-localized gravity 

As in the bigravity models we have two possible ways of realizing a qnasi-localized 
scenario in six dimensions. The hrst model can be built from the a = 0 bigravity 
model, by cutting the space at a point pb on the right of the position of the bounce 
and sewing it to flat space on the right (see Fig.|5.10|). In that case we will have a system 
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P 


Figure 5.10: The cc = 0 quasi-localized model warp factors a{p) and 7(p). The bulk 
region for p > pi, is flat. 


of two positive tension branes and by appropriately tuning the different components 
of the brane tension of the second brane we can hnd solutions where the warp factor 
is constant in the flat bulk region. The hne tunings which achieve this are: 


= — ^tanh 
® 5k 


4A 


^k{pb - po) 


w(2) _ ^w(2) 
^0 ~ 


(5.33) 


From the above formulas it is obvious that if we paste the flat space at the position 
of the bounce, i.e pb = po then the second brane is tensionless. If we relax these hne 
tunings we can hnd solutions where the warp factors in the hat bulk region are non 
trivial, but we will not discuss this possibility here. 

The above construction will obviously have no normalizable zero mode due to the 
inhnite volume of the system and the KK spectrum will be continuous. However, as 
in the hve dimensional case we expect that the low part of the KK tower will have 
a “resonance”-like coupling to matter on the "-I-" brane at p = 0. This will mediate 
normal gravity at intermediate distances but will change the nature of the gravitational 
law at ultralarge scales. 
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P 


Figure 5.11: The a = 2 quasi-localized model warp factors a{p) and 7(p). The bulk 
region for p > p;, is flat. 


The second model of quasi-localized gravity can be built in a similar fashion from 
the a = 2 bigravity model (see Fig. |5.11|) . Again we obtain a system of two positive 
tension branes which will have a constant warp factor in the flat bulk region when the 
following hne tunings are demanded: 


1/i^^ = —^tanh 
ok 


4A 


■^Kpb - po) 


y(2) _ 3^(2) ^ 8^" 

*^0 ~ Q ^9 + 


5p2 y{2) 


(5.34) 


From the above formulas it is obvious that if we paste the flat space at the position 

f2) (2) 

of the bounce, i.e pb = po then the tension Vq vanishes whereas the tension V() 
diverges. This limit is not well behaved since the internal space becomes singular in a 
whole hne. At this limit obviously classical gravity breaks down and one would expect 
that quantum gravity corrections would resolve the singularity. However, one should 
check if these quantum gravity corrections affect the solution in the vicinity of the 
four-brane at p = 0, e.g. by examining the correction to the Einstein action. 

The above construction will again have no normalizable zero mode and the KK 
spectrum will be continuous. A “resonance”-like coupling of the KK states to matter 
on the brane at p = 0 is again expected which will give normal four dimensional 
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gravity at intermediate distances but will change the nature of the gravitational law 
at ultralarge scales. 

3.2 The crystal universe model 

Another obvious generalization of the hve dimensional multigravity models is the one 
of the crystalline brane model. All that one has to do is to paste an infinite array 
of bigravity models. We again take this constructions in two copies. One for a = 0 
(Fig. |5.12| ) and one for a = 2 (Fig. |5.13D . A band structure is again expected as in the 
five dimensional case and the width of the first band will be exponentially smaller than 
the one of the first forbidden zone. As far as the phenomenology is concerned, these 
models will generate normal gravity at intermediate scales as the first band will behave 
as an effective zero mode, whereas at ultralarge scales will have the same change of 
the gravitational law as the one of the corresponding quasi-localized models. 
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Figure 5.13: The a = 2 crystalline model warp factors a(p) and yip). 

4 The weaker energy condition 

In this chapter we have constructed, for the hrst time, flat brane theories which can 
lead to multigravity models and their associated modihcations of gravity at large 
distances without introducing moving negative tension branes. The constructions are 
made possible by going to a six dimensional theory. In hve dimensions with flat branes, 
the presence of a bounce of the warp factor was linked to the violation of the weaker 
energy condition. This is not true, however, in six dimensions as we will show in this 
section. In hve dimensions with the metric: 

+ dp^ (5.35) 

one can readily show that the weaker energy condition requires that: 

A" > 0 (5.36) 

which is violated at the position of the moving negative tension branes. In the case 
that the branes were AdS one could have a bounce (without having negative tension 
branes) and still satisfy the weaker energy condition because the above relation is 

















4 The weaker energy condition 


93 


modified to [19]: 

A" > -2H^e^ (5.37) 

However, such models do not lead to modifications of gravity at large distances and 
moreover the remnant negative cosmological constant is in conflict with current ob¬ 
servations. 

In the six dimensional case with metric: 


ds^ = e + dp^ -|- e (5.38) 

from the weaker energy condition one hnds two inequalities which can be cast into the 
following relation: 


1(5')2 _ Ib" - -A'B' < < --(B'f + B” 

6 3 6 2 


^A!B' + 2{A!f 


(5.39) 


These inequalities hold for the a = 2 models as long as the three-branes sitting on 
the conical singularities have positive tension. It is violated, however, in the a = 0 
models everywhere in the bulk. This shows that in six dimensional flat brane models, 
the presence of bounces of the warp factor is not necessarily linked to whether the 
weaker energy condition is satished or not. It is not yet clear if the violation of the 
weaker energy condition in the a = 0 models is a sign of a possible instability. 



Chapter 6 
Conclusions 


In thesis we have presented a comprehensive account of the multigravity scenario and 
how it can be realized in warped brane-world models in five and six dimensions. It 
is clear that the possibility of having massive gravitons contributing significantly to 
intermediate distance gravitational interactions is not excluded by current observations 
and furthermore it offers phenomenological signatures which are in principle testable. 
Moreover, taking into account that massive gravitons necessarily arise as KK states 
from the dimensional reduction of a higher dimensional theory, these modifications of 
gravity can be a very unusual window to extra dimensions. Multigravity models which 
involve brane configurations with interbrane distances of the order of the Planck length, 
modify gravity not only at distances of the order of the Planck length as usually do 
ordinary KK theories, but also at ultralarge distances. The crossover scale of normal 
to modified gravity is related to the interbrane distances and in that sense is generated 
by the geometry itself. 

The hrst attempts to build such a model faced inconsistency problems. The "H- V” 

bigravity model, the quasi-localized GRS model and the crystalline model, although 
they had promising phenomenology, they all shared the characteristic of moving nega¬ 
tive tension branes. As we explicitly discussed in this thesis the moduli corresponding 
to the fluctuations of these branes are ghost helds and therefore unacceptable. How¬ 
ever, the general characteristics of these models persist in more involved constructions 
that have less or no problems. For example, the appearance of light states whenever we 
have more than one positive branes that localize gravitons is a generic characteristic 
which persists in the AdS brane models or the six dimensional ones. 

The AdS models on the other hand, although they do not have the ghost radion 
problem, do not give some interesting phenomenology. The first reason is the very 
fact that the remnant cosmological constant is of negative sign, in contradiction with 
observations. The second reason is that they do not lead to any observable modifi¬ 
cations of gravity at any observable scale. The modifications of gravity at ultralarge 
scales are hidden behind the AdS horizon, while the modifications at all scales due to 
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the different propagator structure of the massive gravitons are enormously suppressed. 
It is, however, interesting how the AdS brane models escape the van Dam-Veltman- 
Zakharov no-go theorem about the non-decoupling of the additional polarization states 
of the massive graviton. This in addition to the observation of Vainshtein that the 
latter no-go theorem may not be valid even for flat background, suggests that the 
massive graviton proposal does not face any fundamental obstructions. 

Certainly the only theoretically and phenomenologically viable possibility of multi¬ 
gravity at the moment is the one in six dimensional models. Due to the non-trivial 
curvature of the internal space one can have bounces of the warp factor with only 
positive tension flat branes. There is no five dimensional effective regime in these 
models, since otherwise one would face all the problems present in five dimensions. 
The low energy theory is directly four dimensional. Moreover, the replication of the 
multigravity models in two types each, one with conical singularities and one without 
them, gives a wealth of possibilities that offer promising phenomenology. 

There are a lot of open questions regarding multigravity at the moment. The 
most important one is the construction of a cosmological model of multigravity. So 

far the models we considered involved either flat or AdS branes which are not valid 
descriptions of our Universe at large scales where the modihcations of gravity due to 

the massive gravitons are expected to appear. A cosmological generalization of these 
models would probably reveal corrections to the four dimensional Friedman equation 
at cosmological scales. This information is essential if one wishes to put these models 
into observational test. 

Furthermore, as far as model building is concerned, there are a lot of issues to 
be investigated in the six dimensional models. A careful treatment of the angular 
excitations in the case of the non-singular models should be carried out to examine 
the nature of gravity on the four-branes. Furthermore, one should examine the moduli 
of the system, the radions and the dilaton, and calculate their mass. It is quite probable 
that the modulus corresponding to po is massive as in the five dimensional AdS 4 branes 
case and so the system is self-stabilized. Finally, the non-trivial way that locality is 
preserved by these model should be carefully studied. 



Appendix A 

Effective action for 
graviton-radion-dilaton system 


It is convenient to define a new variable z defined by: 

1 dz 

a{y) dy 


(A.l) 


In the coordinates (x, z) the metric ( |3.1CI|) is conformal to a flat perturbed metric 
Gab = Vab + Hab- 


ds^ = m 


dx^dx'" + Gzz dz^ 


= a 


r]MN + Hmn dx^dx^ 


z) + 2e{z) d^dyf 2 {x) + h^^{x, z) 
Hzz = ^2{x,z) 


(A.2) 


where the helds tpi and <^2 are: 


ipi{x, z) = Q{z)fi{x) + B{z)f 2 {x) 
(P 2 ix, z) = q{z)fi{x) + A{z)f 2 ix) 


(A.3) 


Inserting (A-2|) in ( p.l|) and taking into account the equation of motion satished by 
a one gets: 

^3 r 


Seff = / ax Ceff = / ax 2M'^ / dz 


a^cfpQi) + ^ \{d,hf - dzK^ 


Gin — G\ + G2 + G 


12 


+ + Gj^^ 


(A.4) 
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where 


^hip — 




a" I + /2^ (e'a=^) 


dpd^h^'^ 


3 d 


Oh + -— {a^a'ip2 - a^T^i) h 


(A.5) 


-^^1 = df,fid^fi (Q^ + Qg) + + Sa'^aq^ - Ga^aQ'q^ (A.6) 


^^2 = dj2d^f2 ( + ^a^AB - Sa^B'd + Sa^a'e'A 


+ fl {^a^B'‘^ + 3a'^aA2 - Qa^aAB'^ 


(A.7) 


Ci 2 — d^fid^f 2 ( —aPBq + —a^AQ + Sa^BQ — Sa^e'Q' + ?>a^a'e'q 


/1/2 


Qa?B'Q' + Qa'^aAq — Qa^a' {AQ' + gi?') 


(A.8) 


and C^p\(h) is the four dimensional Pauli-Fierz Lagrangian for h without the mass 
term: 


^%Ch) = \d,K» - \djidih<^^ + -^djid^h (A.9) 

All indices in the above expressions are raised and lowered by rj^y. Differentiation 
with respect of z is denoted with a prime. The absence of mixing terms in Dg// yields 
the following constraints: 


A ( \ 

^ ^ = — , 
a 

4- (So'O + A-4- {ah') = 0 

dz ^ ^ dz ^ ^ 

(A. 10) 

Q{z) = ca~'^, 

q{z) = —2ca“^, c is a constant 

(A.ll) 

J dz aA{z) = 0 


(A, 12) 










Appendix A: Effective action for graviton-radion-dilaton system 


Equations ( |A.1(J|) - (|A.12|) give: 


£,ff = 2M^ / dz-f^a^E^^lCh) + j [{djif - 




(A.13) 


In particular the effective Lagrangian C-scai for the dilaton /i and the radion /2 is: 


^Scal /Ci/iD/i + ZCa/aD/a 

/Cl = 2M^ -c^ f a~‘^ dy 
2 J-L 


with: 


4- (Ba^) + 2a 
dy 



da d 
dy dy 

dB 

dy 


{a^dyc) 
-- 0 


0 


(A.15) 


As a result only dye{0), dye{r) and dye{L) enter the radion effective action. 









Appendix B 

Ricci tensor and scalar for the 
dilaton perturbation 


In this Appendix we list the Ricci tensor components, the Ricci scalar and the action 
obtained by the metric: 

^'^^^^a{yYg^y{x)dx^dx'^ + (^ + (B-1) 

The spacetime components of the hve dimensional Ricci tensor are: 


dig gi/ 


9ai/ —2r-i du 


a ^7 A _ 1 1 - g ^7 '\ -4 

2 (l + a-27)^’''^ 2 VI+ 0-27/“ 


—4 

a 7 


(1 + a-27) 




3 + 4a 27 
1 + a~‘^'y 


9^lu^ 


—a ^7 


aa 


(1 + 0-27) 


(B.2) 


and the (55) component: 


R55 = e“''^(l + a-7a-7^7’^-e“''^(l + a-7a-^n7 
—4(1 + a“^7)a“^a" — 4(1 + a“^7)a“^(a77 


(B.3) 


Finally the Ricci scalar is: 


R 


^^a-2 


Rg + e“''^ 


A + 3a 277 
V 1+ a-27 j 


‘□7 + 


A — 3a 277 

V 1 + a-27 J 


a ®7,^7’ 


M 


a ^a" 


(1 + 


12a 2(a')^ + 20a '^(a')^7 

(1 + 


(B.4) 


In the above expressions the indices are raised and lowered by 
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Appendix B: Ricci tensor and scalar for the dilaton perturbation 


The action (|3.48|) then becomes: 


S= d'^xdyy/^ < 2M^ 


e “ "V(l + a '^'^)Rg + e “ (1 + ^7 + V ) 


+e“2“ '^(-8a"a^ - 12(a')^a^ - 20(a')^7) 


1 


2(1 + 0 - 27 ) 


e-2“ "7$/. _ e-2a ^ 




Ai('t)i(9 - k) 


(B.5) 







Appendix C 

Equations of motion in six 
dimensions 


In this Appendix we provide the Einstein equations for the class of models that we 
considered in chapter 5 where the metric is: 

= a{p)ri^ydx^dx‘' + dp^ + 'y{p)d6‘^ (C.l) 

Ignoring the three-branes located on the conical singularities, the Einstein equa¬ 
tions are generally written as: 

+ rS) (c.2) 


where is the bulk energy momentum tensor of the general form ( |5.12|) and 
the four-brane energy momentum tensor of the form (5.13|). 


In the absence of a four dimensional cosmological constant, the {6, 9) component 
of the above equation is: 



a 




Ag 

4M4 


AM^ 


- Pi) 


(C.3) 


The (p, p) component is: 


3 /^ ^ Ap 

2 \ a J cry 4M^ 

Finally, the (p, u) component is: 


(C.4) 


3 a" 3 cr'y' 

o-^ 7- 

2 cr 4 cry 






S{p- Pi) 


2 


4M4 4M4 


(C.5) 
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Appendix C: Equations of motion in six dimensions 


These equations may be compared with the ones of the hve dimensional case where 
the metric: 

ds^ = a{p)rjij_udx^dx'^ + dp^ (C.6) 

gives rise to the (p, p) component: 


2\aJ AM^ 


(C.7) 


and the (p, z/) component: 


3 a" 

2V 


A 

4ii^ 


1/* 

4ii^ 


<5(p - pi) 


(C.8) 


The extra freedom that we have in the six dimensional case is apparent. 
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